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In this Issue: 

I Our cover photograph shows a HP 77020 A Ultrasound Imaging System being used to 

display a picture of a patienfs beating heart. Ultrasound imaging— bouncing a series of 

high-frequency sound pulses off the organs and other structures inside the body — lets 

j^ cardiologists^ obstetricians, and radiologists see the body's inner workings without surgery 

I iT ' ^^- ^^^ without any health risk to the patient. Heart valves and chambers can be seen in motion. 

-'^^^P ^ developing fetus can be obsen/ed and measured without risk. The HP system, we are 

"1^^'^^ told, is particularly easy to use and produces an especially high-quality image. There's no 

J!r"'^!^ I question that an enormous amount of engineering talent and effort went into its development. 

It's designers gave us so much interesting information that it will take nearly two full issues to publish it alL 

The articles in this issue discuss the system and its use, describe the design of Its novel transducer (the part 

that touches the patient), and relate what was done to ensure image quality. On page 13, a cardiologist gives 

us the user's point of view. In December, we'll complete the story with articles on the hardware and software 

design, 
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Ultrasound Imaging: An Overview 

By using a beam of ultrasound, it Is possible to look at 
organs and other structures inside the human body without 
breaking the skin, 

by H. Edward Karrer and Arthur M. Dickey 



WHAT IS AN ULTRASOUND MAGE? An image 
is a reproduction or imitatioD of the form of 
something. A major thrust of modern technolog}^ 
in the past century has been to enhance or complement 
the capability of the human sensory system with new im- 
aging modalities. Examples that come to mind include the 
optical telescope, radar, sonar, mfrared, and X-ray imaging 
systems. 

The imaging process consists of two parts. First is the 
illumination of the objects to be imaged by some radiation. 
This radiation interacts with the object by means of reflec- 
tion, absorption, or scattering. The second step is a recon- 
struction of the energy received h:om these interactions to 
form an image of the object. 

The primary goal of medical imaging systems is to '*see" 
inside the human body noninvasively. Our eyes are not 
useful for this since biological tissue is largely opaque to 
visible iigbt. Fortunately, other forms of radiation like X- 
rays, some nuclear particles, and ultrasound readily pass 
through and interact with tissue. 

Ultrasound brings several useful capabilities to medical 
imaging. An ultrasound image represents the mechanical 
properties of the tissue [i.e., parameters such as density 
and elasticity]. Most people can immediately recognize 
common anatomical structures in an ultrasound image 
since the organ boundaries and fluid-to-tissue interfaces 
are easily discerned. This mechanical information comple- 
ments the results of other imaging techniques, such as X-ray 
imaging, which displays X-ray absorption. 

The ultrasound imaging process can be done in real time. 
This means that the viewer can follow rapidly moving 
structures such as the heart without motion distortion. 
Also, it allows the operator to move the small handheld 
transducer probe over the surface of the patient's body to 
select the proper view in an interactive way. 

Ultrasound appears to be one of the safest diagnostic 
imaging techniques. It does not use ionizing radiation hke 



X-ray and hence is used routinely for fetal and obstetrical 
imaging, There have been many studies done on the possi- 
ble biological side effects of nltrasound on tissue and none 
has shown any harmful effects during clinical examination. 
Areas where ultrasound imaging is particularly useful 
include cardiac structures, the vascular system, the fetus 
and uterus, abdominal organs such as the liver, kidneys, 
and gall bladder, and the eye. However, ultrasound cannot 
be used to image all portions of the body. Air pockets, for 
instance, are excellent reflectors of nltrasound and limit 
the penetration of sound into the lungs and bowels. Bone 
highly attenuates ultrasound, which makes imaging the 
adult brain through the skull difficult. 

Physical Principles 

Ultrasound generally refers to sound that has a frequency 
or pitch above the range of human hearing. Sound waves 
in human tissue are compressional (longitudinal) waves. 
These waves consist of repetitive or periodic regions of 
local compression and rarefaction which propagate 
through a medium at some velocity v. The spacing between 
successive regions of compression is called the acoustic 
wavelength X. These parameters are related by the expres- 
sion v = fx, where v is the velocity in meters/second, f is 
the frequency in hertz, and X is the wavelength in meters. 
Because it is not possible to resolve objects smaller than a 
wavelength, it is desirable to image at the highest possible 
frequency (smallest wavelength] to get the best resolution* 

The most common mode of ultrasound imaging is the 
reflection or pulse echo mode. The analogy to radar is valid. 
A short burst of acoustic energy is directed into the body 
and the reflected energy is received at a later time. Different 
tissues reflect acoustic energy depending on their charac- 
teristic acoustic impedance. This impedance is defined by 
the expression Z = pv, where Z is the acoustic impedance 
in rayls, p is the density of the tissue in kilograras/cnbic 
meter, and v is the velocity in meters/second. 
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Fig. 1 . Phased -array operatfon of 

an ultrasound imaging syfitem. By 
varying the delay Ti for each trans- 
ducer element electronically , the 
direction of the outgoing ultra- 
sound pulse can be steered and 
the array can be dynamically 
focused on the returning echoes. 
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As the ultrasound pulse passes from one tissue type to 
anothert a portion of it is reflected in a way analogous to 
reflection from a discontinuity on an electric transmission 
linen For instance, when a pulse crosses the boundary be- 
tween muscle tissue (Z = 1.7X10^ ray is} and blood 
(Z — 1.6x10^ ray is), about 0.1% of the acoustic energy is 
reflected. This echo is used to create the image of the blood- 
muscle interface. The rest of the pulse propagates across 
the interface and continues into the body to image deeper 
structures. 

The ultrasound pulse is transmitted and received with 
a piezoelectric transducer. The piezoelectric material is a 
poly crystalline ceramic which converts an electrical signal 
into acoustic energy and vice versa. This transducer is the 
source and eye of the imaging system. It is usually placed 
directly on the sldn surface. 

As the ultrasound pulse propagates through the body, 
its acoustic intensity is attenuated with propagation dis- 
tance. The acoustic attenuation also increases with fre- 
quency. A typical attenuation coefficient for soft tissue is 
1 dB/cm/MHz. When imaging a structure that is 10 cm deep 
in the body at 3.5 MHz, the round-trip attenuation is 70 
dB. This is an easily detectable signal when using clinically 
acceptable transmitter powers. Note that at 10 MHz the 
round-trip attenuation would be 200 dB. This points out 
a basic tradeoff in ultrasound imaging; as the frequency is 
increased to resolve finer structures, only shallow struc- 
tures can be imaged. Typically. 2,5 to 3,5 MHz is used for 
deep abdominal imaging, 5 MHz is used for organs near 
the surface and in pediatric work, and 5 to 15 MHz is 
reserved for "small parts'' imaging such as the thyroid, 
testicles, peripheral vascular vessels, and the eye. 

Other structures can also interact with ultrasound, giving 
different tissues different textures. Objects much smaller 
than a wavelength cause Rayleigh scattering. This diffuse 
scattering depends on the fourth power of the frequency 
and the ultrasound pulse tends to scatter uniformly in all 
directions. Images made at higher frequencies tend to ap- 
pear more '^filled in" because of this effect. If an object is 
much larger than a wavelength, it becomes a specular re- 
flector and acts much like a mirror. 

Ultrasound Imaging Systems 

B or brightness mode scanning is the most common type 
of ultrasound imaging, A B scan is a view of a cross-sec- 
tional slice through the object. A narrow pencil beam of 
ultrasound is swept through a sector to define the scan 
plane. The beam is formed from bursts of ultrasound. The 
repetition rate of ultrasound pulse generation is selected 
so that the transmitted pulse has time to travel to the 



deepest target and back again before the next pulse is 
launched. The pulse is assumed to travel in a straight line 
at a constant velocity, a good assumption in soft tissue. As 
the pulse propagates into the body along any scan line, 
echoes are generated which travel back to the receiver. These 
echoes vary in intensity acording to the type of tissue or 
body structure causing them. This data is presented on a 
cathode-ray-tube (CRT) display in which the brightness 
(hence the term B mode) depends on the echo strength. 
This image can be formed in real time since it takes only 
260 piS for ultrasound to travel 20 cm into the body and 
return. If the display frame rate is 30 frames per second, 
125 lines can be displayed in the image. 

Ultrasound imaging systems differ in the way that the 
sector scan is implemented, A single transducer can be 
mechanically rocked or pivoted to create a sector scan. The 
phased-array approach uses a stationary array of many 
small transducers which are electronically controlled to 
steer and focus the beam. This concept is shown in Fig, 1. 

In effect, the phased array acts like an acoustic lens with 
electronically variable focal length. Assume that a region 
of tissue at location O has reflected a small amount of an 
incident short pulse of ultrasound. The echo reaches each 
transducer element [TJ along the array at a different time 
depending on the position of O. Variable time delay ele- 
ments (tj) in series with each transducer element are 
rapidly varied under computer control to bring the electri- 
cal pulses generated by the echoes received at each trans- 
ducer element into coincidence at the summing junction. 
For instance, if the object is on the axis of symmetry of 
the transducer array, but is far away, the returning echo 
wavefronts are essentially planar and the time delays are 
set equah If the object is off-axis at some angle B in the far 
field, the time delays must be linearly varied across the 
array aperture for coherent addition at the summing junc- 
tion. When the object is close to the array, as shov^m in Fig. 
1, the approaching wavefronts are more curved and the 
delay distribution across the array is more nearly parabolic- 
Changing the individual transducer delays to accommodate 
the curvature of the wavefronts is called dynamic focusing 
and is necessary to achieve maximum resolution in the 
near field of the array. 

In a similiar manner, the outgoing ultrasound pulse can 
be directed to a desired location. By varying the time at 
which each element of the transducer array is excited by 
the transmitter, the direction in which the ultrasound 
pulses generated by each element coherently add to form 
one pulse can be changed from one side of the axis to the 
other. This is called beam steering. 
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Fig, 2. Block diagram of a basic 
ultrasound imaging system using 
a phased-array transducer. 
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History of HP's Ultrasoimd System 
by John T, Hart 

RSfD Manager. Andover Division 



n j hope 10 convey me S9n% of ac^mv^rwit^ and sorrsimi^ cSpjiijsionrRerif that w^ 

Wien Ed Kafrer's groue ^ HP Labofaiones (MPk) ipc*i up me c^iange sf buffing 
a phasec^^iTay system a numt^f of ye^s ago. ntrm expected mat the HP 77Q20A 
would be tie msuA. Dave Perczoii, noMr genera^ managef of the Andoffir DIvffiian. 
pro^rldedtheinitiflllfiaisonformedMekin. ^^Dr RictiardRipprrQmSiaifordUnhi^orsHy 
HoGfMtal provided the ciinicaf qui<iai>oe K^any Layers of Itie onion t'cad ^ De r^noved 
balore we mei our perfofmarce gaai The ffr si system was impressed, and occuphed 
^lalf the area of a large fOOfti, out tnt reaction to Hie images was genuine ewcstement 
Orme TV screefi was a picture Gf John Larsoo^s twaiing (i^rt 7Ta.( was the umpetus 
to lorm a seonon at the Andover D'VfSion to develop the techftoiogy arvd txjiSd a new 
product iine- 

WWIe the first images provided a very real (ikenesfe, they were o¥e^ayed with "ghost ' 
imege*. This artifact was caused by gratfng tobes, in turn caused by sparse tfansdiicer 
aiffiy apaqing. sirvce there *efe onty 16 elerrients To eliminate this, it was necessajy 
to have 64 elements in the iransducer, spaced one haJf wave^engiti apart This was 
^scouragirtg because i! was just too costty to have that many receiver channels^ each 
with Its own deJay i^ne Fortunately, however. Sam Maalalk had dessgned a mixing 
capability into the first system whrch wou^ct al'ow srna<l time delays to be appfoxhmaled 
by phase shifta ir two months, it was tried ar>d found to work It was possible to build 
a S4'e»emenE system el reasonable cost and eliminate the graiiria: lobes oausmg the 
ghost umages 

h WH& at this point that the promise of a high-perfcrniahce cardiology system at 
subsEantiaily lOwer cost Pecame apparent Ther^ were a few real-time ultrasound 
systems in ttie marketplace at thai time which clearly demonstrated the advantages 
of leal-lime uJirasound. but they were expensive and the picture guasity was far from 
optinoal. With lols of determination and the promise of a reai: advance, we undertook 
a jotnt HPUAncTowBr Division effod to build a new system. 

The iniiiar Andower organization was small* Arfriur Dickey was p^ojact manager for 
the scanner, transducer and scan cdr^verier Larry Banks wa& responsible for system 
Integration and control As the profect grew in S42e, Jim Feamslde took on the task of 
managing the transducer dsveiopment, Ray O'Connell ihe display and scan converter, 
aid Tony VaJlance the system software I was the section manager 

Sam Maslak worthed wrth the scanner team to design Ihe beam sleering electronics 
This was 10 be driven by HP'g custom te-bit silicon -wi-sapph ire microprocessor instead 
ot tha orjgmal mjmcomputet HPL had lo carry the tjall m developing the 64-83emen! 
iranaduce^ as Andover had neither the facilities nor the expertise to do it The highest 
poorcty was to develop the scanner and iransducer. Because that was wt^re the basic 
tmprsvement m image quality would be denvetd AJso. our theories had lo t>e tested 
in praclioe After considerable work, there was iniftial success, mixed wath more disap- 
polnimenL The images showed ihaE the grating lobes were eirmineted, hui ihe images 
were lilled with clutter and had lower resolution than expected 

The solution required extensive reshaping of the acoustPC bearr eiiminatkm of spurn 
ouB transducer modes, and careful controf of (he frequency ^v^ '^unn Hugh Larsen 
and Rick Pering were able to modal the design parameters ario -..alGuiate the new 
design constraints An Dickey and Rori GatifSte rose to the ctislliBnge, squeezing out 
aignllicant partormence improvements ^n the scanner, while flob Mc Knight was given 
tha ditticuii assignment of optimizing the video processing chain Jim Feamside and 
John Larson provided some innovative soluljons to our transducer problems, white 
Dave MiJfer contributed m line complex area of acoustic lens design 

Tha image quality got another boost from Hugh Larsen and Steve Leavitt wtio 
invented Ihe R-Theta scan conver^Jon algonthm Barry Hunt joined in This effort at a 
later time and made significant conifibutions. This sctieme nnakes it possible lo acfiieve 
artifacf-lree reoonstnuction of the ultrasound image ^n television video output. The 
hfyquisi criterion as 11 relates to the aooustic fine density Indicated that this could be 
adiieved, but up until this time all reaJ-time sector ultrasound images were plagued 
with a stappy or lagged appearance because of the neces:&ariiy simple approximations 
made in cafcuiatmg a picture "on the lly." The 77020A system demonstrates a display 
!n whk:h such discontinuities are atmosi irnpgrceptibie 



In inraf^ vAth trsssfiod 10 ifnpfQweiTaoe equality. Ptei^e Rhtoa^ 
^stpi. and L^ry SanKs <*i^m a^if^iJnQ the p^^taf^ design and t%jf^n tnierfa^ 
to Giinical trials Hm^. ih& 6D-Hr CRT dbpkary amd the pehfiiera^ device such as 6m 
vJOeo casaetse recofder and Ihe stn^pc^iaii reooroar prmefmfi e pas^aMai^ k -n potimn 
dmli&i^. becsuse fbe Images thai the is& sees rnost come through all o! these 
devices it turTwd out that rione <if these p>robiems was e^y to solve, because tt« 
eye is arr^zmoJy s^isitive to &tein the sl^hteet irtaoz^dcy m gr^ scale of position 
John Duh^ at HPL had disooviefed some d ijie gray-scaie rrvapping pTobiems m the 
earlier HR system and had dorie some pretimrnaiy wo^tc to quanti^ly these cornpjex 
problems Al Tykutsicy, assisted by Jan Acoettura, worked out the problems in the 
60- Hj deptay Steve LeavJtt Jim Mnece, ana later Afwyn D'Sa aoived the problems 
m the ODgiiai v<oeo miertaces. Jim Corn ad and fater Rjch Jundanian worked out the 
interfaces mvotved in putting mages on the sm^s^hart recorder 

The human intertace also received a great deal of attentiori dunng the Clinical tnats 
Within the 77O30A, there are muHnpie l&-bit microprocessors, each operating its own 
subsystem I] is imponant in a clinical environment that the operation be obvious ar>d 
naluraf. and that the system recover from errors easily This presented a significant 
challenge to the software design teams Joe Luszcz and BiN Koppas made significant 
conthbutions to the overall sottwaje system design and the impiementation of an 
effectcve and reliable reaf-time operating system Definition and «mplementatJon of Ihe 
scan converter architecture was done by Dave Hempstead, while delmrtion and ittv 
plementation ot the scanner software was done by Al Langguth and Bob Kunz. The 
high-speed communication of data between these subsystems requfred careful defi- 
rtition ol archiiecture and communication protocols and was an excel leni ejsampte of 
the cooperative team eflort that characterized the proiect 

Testing and Servicing ftie 770^A System were equally complex and challenging 
because so many new technologies were being mtfoduced Simultaneously. Radhu 
Sasu, Tom James, and later Tony VaSlance and A I Langguth each made feat contribu- 
tions lo the tesi plan and m managing us implementation Tony Fisher designed an 
innovative auiomatiic transducer test system Idr beam plots and individual aiomant 
characterization, permitting Ihe maniLifacSure of consistently high -quality tfandsucers. 
Gary Sea vey and Tom S;jabo worked out many of the probiems involved in characieriz- 
ing the transmitted acoustic power This ts especially important, since the "object" 
under examination is a human bemg, 

Continued work has further Improved system performanoa stnoe its initial Introduc- 
tion Several 3.5'MHz and 5 0-MHz transducers, which significantly improve image 
quality, have tieen added lo iha product line Quarter- wave matching technotogy has 
been developed to increase transducer sensiiivify. permctfing these probes to be used 
on most patients Jim Fearnside Dave f^iHer. and Jerry Leach p'ayed tha key roies in 
developing these transducers Tom S^atxi developed a more accurate model for 
transducer eSements which speeded up the devesopmont ol these ttansduceis and 
gave us a toot that accurately predicts ttie iperformance o( new transducer configura- 
flons 

Hon Gatzke's scanner group, notably Ride Shyder, Syd Karp. and Jim Conrad, save 
image quality another boost by Increasing tha number of transmitters to 64 and reshap- 
ing the beam and the gray scale A new image processing team, headed by Paul 
Magnin, has tieen chartered lo find even more ways to improve the very tmportant 
parameter of image quality 

AJong with higher performance, a broader range of penpherals and system enhance- 
ments such as Rachel Kinjchi's quaniiEJative cardiac and obstetric analysis packages 
have provided a high-perfonnancB. artractirve, mobile, and easy-to-use system 

There were many others without whom the project could not have been oompieted. 
In all, itiere v>if©re well over tOO people actively involved, irtcluding the printed circuit 
design and fa&riciatFon teams and the manufaciuhng engineers wfw were er^aged in 
(he desi-gn and development of Ihe processes and the product 

The HPt and Andover Division development teams have designed the hlgh^t-per- 
formance cardiology ultrasound system available today, and perhaps more important, 
have put into place a solid base of technology and application understanding on which 
we can extend the use ot diagnostic ultrasound into new areas of he^fh care 



In practice T the transmit and receive functions are sepa- 
rate. Consitier the events necessary to generate one line of 
the image. First, a set of ultrasound pulses is launched 
so as to provide a coherent wavefront along a given scan 
direction. This beam is focused at midrange. As the trans- 
mitted wavefront propagates into the body. It is scattered 
by objects along the scan line and echoes start returning 
to the array. In the receive mode, the array continuously 



changes its focus as echoes are received from successively 
deeper regions. This dynamic focusing process allows an 
array with an aperture of 2 cm to have a varying depth of 
focus fro 01 2 to 20 cm. 

Resolution is the ability to separate small objects in the 
scan plane visually. It has three components: range reso- 
lution (along a scan line), azimuth resolution (perpendicu* 
lar to a stjan line within the plane of the sector scan), and 
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elevation resoJution (the thickness of the sector scan slice). 
The range resolution is determmed by the ultrasound pulse 
length. A shorter pulse gives higher resolution- The pulse 
lengtli is primarily determined by the system bandwidth. 
Two targets 0.75 mm apart In range can be resolved with 
a 1-^s ultrasound pulse (3 cycles at 3.5 MHz). 

The resolution in the azimuth direction depends on the 
array aperture (length of the array) and the acoustic 
wavelength X. The principles are the same as for an optical 
telescope. High azimuth resolution is achieved with a large 
aperture and a short wavelength. Dynamic focusing op- 
timizes the azimuth resolution of the imaging system. 

A block diagram for a general ultrasound imaging system 
is shown in Fig, 2, The transmitter excites the transducer 
eiements with short electric pulses so that a burst of ul- 
trasound is generated. The returning echoes are applied to 
a variable-gain stage in the receiver called a time gain com- 
pensation (TGC) amplifier. This amplifier increases its gain 
witli time, thus compensating for tissue attenuation as 
echoes come from deeper regions of the body. In a phased- 



array system, each element of the transducer has its own 
TGC amplifier The beam former combines the outputs of 
the individual receiver channels by using variable time 
delay and phase adjustment to bring the received signals 
into coincidence and hence bring an object i nto focus . Since 
the received echoes have a very wide dynamic ranget the 
signal processing stage uses signal level compression. For 
instance, when imaging the heart, the strong echo from the 
posterior wall of the left ventricle can be 50 dB greater than 
the diffuse echoes from the endocardial tissue. This data 
must be compressed to be compatible with the 25-dB 
dynamic range of normal CRT displays. The compressed 
signal is then converted from analog to digital form so that 
an entire image frame can be stored in the memory of the 
digital scan converter. The digital scan converter changes 
the scan format from a sector display to a conventional 
television raster scan. It also permits digital postprocessing 
of the images. These various activities — the transmitter, 
TGC amplifier, beam former, scan converter, and display — 
are coordinated under microprocessor controL 



An Ultrasound Imaging System 

This instrument views t/ie internal organs and tissues of the 
human body in real time by di rooting a beam of short 
ultrasound pulses into the body and then receiving and 
prooesslng the acoustio echoes to form a displayed image. 

by Lawrence W. Banks 



THE HP 7 7020 A ULTRASOUND IMAGING SYSTEM 
(Fig. 1) is a real-time phased-array imaging system 
offering a 90* sector image constructed with a polar- 
to -rectangular conversion algorithm that minimizes ar- 
tifacts in the image. The operator can select either a sector- 
only display or a small-sector display with one or two 
M-mode traces. The system provides the operator with a 
variety of capabilities which include stop action, hard-copy 
reproduction, and quantitative analysis. 

Product Goals 

Developing a system of this size and complexity requires 
a clearly defined list of goals. From several marketing 
studies came these five: real-tim:e imaging, phased-array 
architecture, high image quality » mobility, and ease of use. 

Without the ability to image and display in real time (30 
H^)p the system would not address the needs of the car- 
diologist to visualize the heart, especially Its valves. Al- 
though not required in radiology applications , real-time 
scanning increases the speed of the examination and leads 
to better and more cost-effective diagnosis. 

The second goal was a phased-array electronic architec* 



ture. Of the three basic real-time imaging technologies, 
mechanical, linear, and phased-array, the last is the most 
complex, but allows the most versatility and combines the 
advantages of the other two. Phased-array systems use a 
small and reliable transducer that can transmit and receive 
acoustic energ}^ through the space between the ribs. This 
gives a sector format vrith a large field of view, allowing 
the full adult heart to be imaged. 

The most elusive and hardest to define goal was image 
quality. Although easy to quantify from engineering con- 
cepts of acoustic beamwidth, signal-to-noise ratio, pulse 
length, etc., defining what a doctor "sees" is far more dif- 
ficult. In fact^ different types of doctors such as car- 
diologists and radiologists view images differently. 

In 1976, state-of-the-art ultrasound systems were either 
early attempts at real-time systems wfith limited image qual- 
ity or brightness (B) scanners. The images from these sys- 
tems looked very different and yet together satisfied the 
needs of doctors at that time. Our goal w^as to exceed the 
image quality of the single-frame B scanners and improve 
the features of the real-time systems. 

Most systems available during this period were large and 
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F\q. 1. HP Mode/ 77020A Uttrasound Imaging System pro- 
vides both M-mode and two-dimensionaf sector scans of 
human tissue for cardiac, obstetric, and abdominal emmtna- 
lions. The system is mounted on a mobile cart in a configura- 
tion for rapid and easy use. 

difficult to move throughout the hospital. To allow use in 
intensive care and coronary care units, the 7702QA had to 
be small and mobile^ the fourth goal. 

An easy-to-use system is also important in medical appli- 
cations. This became the fifth major goal. Many systems 
were '*stacks of boxes" with controls scattered everywhere. 
The 77020A centralizes all controls on panels surrounding 
the display, allowing the diagnostician to concentrate on 
the diagnostic information in the image rather than on these 
controls. 

To guarantee that the system would be easy to use, all 
hard-copy, display » and storage modules have matching 
gray-scale images. This requires proper gamma correction 
for each medium. Although a simple concept^ many copies 
can he wasted trying to get a good hard copy of the video 
image on the monitor. Slight differences in gray-scale 
linearity can subtly but significantly change the appearance 
of the image. 

To guarantee friendliness, the system was designed from 
the outside in. The human interfaces, including display 
formats and control panels, peripheral locations, size, and 
weight were all determined based upon human rather than 
technical requirements. Although this made the engineer- 
ing requirements more difficult, it helped solve the doctor's 
problem. 



System Controls 

The 7 7 020 A System is controlled by the operator through 
three control panels. The vertical panel to the right of the 
display contains the controls needed for normal imaging 
and peripheral control. The sector imaging controls include 
the major controls for generating the ultrasound image. The 
TGC (time gain compensation) controls used to compensate 
for the acoustic attenuation of the human body allow the 
operator to optimize the image by adjusting the gain of the 
receivers for different depths. 

The TRANSMrt and COMPRESSion knobs control the trans- 
mitted power of the system and the video compression of 
the image. The DEPTH selector adjusts the maximum imag- 
ing depth to 4. 8, 12. 16. 20, or 24 cm. w^hile the SPEED 
selector alloiivs an xM-mode trace speed of 10. 25, 50^ or 
100 mm/s on both the display and the stripchart recorder. 
The final imaging control is the FREEZE button used to 
store the current image in the digital scan converter's mem- 
or>'. 

M-mode imaging along a selected sector display line is 
controlled by the group of controls in the lower center of 
the vertical panel. The paddle switches in the center of the 
M-MODE area select the particular acoustic line. Two lines 
can be displayed at once (Ml and M2). The center paddle 
switch allows any 8-cm segment of the two lines to be 
selected as an M-mode "zoom/' The M-mode image format 
is selected by two pushbuttons which toggle the M-mode 
traces on and off. 

The pushbuttons at the top of the vertical panel control 
the various peripherals. The CAMERA button takes a picture 
on the video camera. The FRAME button freezes an image 
and then sends the frozen image to the stripchart recorder. 
The TAPE button enables the video cassette recorder's re- 
cord function- The CHART button runs the stripchart re- 
corder for normal operation in M-mode- The CHART, TAPE^ 
and FRAME functions are repeated on afootswitch, allowing 
hands-off control of these functions. 

The setup controls hidden under a door (shown open in 
Fig. 1] at the bottom of the panel are used to ciistomize the 
system for a particular image. The two SECTOR switches 
set up the system for the type of organ being imaged. The 
VIDEO switch allows either white on black (W/B] or black 
on w^hite (BAV) imaging. The PRE-PROCessing switch allows 
a choice of four different analog video processing options. 
The POST-PROcessing switch selects four different digital 
gray-scale processing options. The TEST switch and five- 
digit hexadecimal display are used for the internal resident 
test software. 

The horizontal panel below the display contains all the 
necessary controls to annotate the image during real-time 
scanning, play back the image from a video cassette re- 
corder (VCR)» trigger the image from an ECG waveform, 
and analyze the imagen The alphanumeric keyboard is used 
to enter patient ID and other data onto the image- 

The controls available on the horizontal panel to config- 
ure the system or analyze images are generally used when 
the patient is not being scanned. The left part of the panel 
includes a joystick to allow the operator to measure lengths, 
areas, and volumes of two-dimensional Images, and 
heights, times, and slopes of M-mode traces. The ▼ and A 
buttons select which end of the cursor is enabled for move- 
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Quantitative Analysis for Ultrasound Imaging 

by Rachel M. Kinfcki 



ft fS often necessary to supplement the quafitative anatomical 
and functiorial information obtained during an ultrasound exami- 
nation with quantitative measuretrients and derived calculations, 
Since cardrologfsts and obstetricians, in particular, rely heavily 
on quantitative intormation obtained from ultrasound images, two 
analysis packages were added to the HP 77020A Ultrasound 
Imaging System to meet their needs, These packages provide 
measurement capabiiity via the system calipers to measure de- 
sired parameters accurately, perform standard calculations 
based on these measurements, ailow for input of patient history 
data and comments, and generate a hard-copy report that sum- 
marizes the data. 

Cardiac Analysis 

The ability to assess ventricular function is one of the most 
important applications of ultrasound in cardiology. To distinguish 
abnormal from normal heart functioning, the left ventricle, in par- 
ticular; is studied thoroughly using measurements and derived 
calculations. Linear dimensions such as chamber diameters, wall 
thicknesses, and valve separations are determined as well as 
cross-sectional chamber areas, Derived quantities, such as left 
ventncuEar volume, stroke volume, and cardiac output, are calcu- 
lated from these measurements. 

The heart's left ventricle pumps blood returning from the lungs 
to the rest of the body. Changes in teft ventncular blood volume 
over the heart cycle are an indication to the cardiologist of the 
ventricle's functioning ability. M-mode and two-dimensional im- 
aging techniques provide one- and two-dimensional measure- 
ments respectively, but volume es a three-dimensional quantity. 
Much research has been done on accurately computing left ven- 
tncular volume from these measurements and this has led to 
modeling the left ventricle as a geometric figure or combination 
of figures. For example, the bulfet formula assumes the ventricle 
is shaped tike a combination of a hem f sphere and a cylinder 
and the modified Simpson's rule formula assumes the shape is 
a combination of a cyiinder, a truncated cone, and a cone. The 
cardiac analysis package offers the operator a choice of sly^ of 
the most accepted techniques for calculating left ventricular vol- 
ume.. 




Obstetric Analysis 

The obstetrician uses ultrasound as a tool to determine fetal 
viability and to monitor fetal development. The gestational sac 
can be identified with ultrasound by the end of the fifth menslrual 
week (fetal gestational age is most often computed from the first 
day o! the mother's last menstrual period). By the seventh 
menstruaJ week, the fetus can be recognized and a pulsating 
heart can be seen. As the fetus develops, ultrasound reveals 
fetal anatomy, placenta development, and uterine structure, 

Much research has been done to correlate fetai anatomical 
dimensions with fetal age and weight. Fetal age determination 
is probably the most important information available to the 
obstretrrclan since it is used to detect normal versus abnormal 
grov\^h. One of the most accurate ultrasound techniques for de- 
termining fetal age is to measure the fetal crown- rump length in 
the first trimester of pregnancy when bioi'ogicai variation between 
fetuses is usually small. The crown-njmp length is the longest 
length of the fetus excluding the limbs. Several researchers have 
measured this length when the gestational age was known and 
have established tables or formulas that can be used to verify 
age from this measurement In the second and third trimesters, 
the crown -rump length becomes too large to be completely seen, 
so measurements such as the fetal bi parietal (head) diameter 
and femur length are used to determine fetal age. These mea- 
surements, together with abdominal measurements, can be used 
to estimate fetal weight and to investigate the possibility of intra- 
uterine growth retardation which may indicate that the fetus is 
not receiving proper nourishment. 

Analysis Packages 

Two phiiosophies of how analysis packages should work were 
discussed when designing the on-line analysis packages for the 
7 7 020 A. The first approach can be called the walk -through 
method; the operator is prompted for ail measurements in a 
predetermined sequence and must either make each measure- 
ment or skip over it before continuing with the next. The second 
philosophy is to allow free-form input of measurement data. This 
method allows the operator lo input measurements in any order. 
It is the operator's responsibility to enter all of the necessary 



Fig. 1. QbstBtfic analysis package measurement menu with 
tela! biparietBi diameter measurement stored ^ 
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Fig, 2, Sample edit report for obstetric anaiysis package. 
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measurements during the course of tfie study The HP anaiysis 
packages are implemented using the philosophy of free-form 
measurement input. Although the walk-through method might be 
better for ine^penenced operators, rt can be very fmstrating to 
the experiericed operator wtx) prefers to make measurements 
m a logical (perhaps anatomical) order that differs from the 
package measurement order This trustratioo is particularly 
rvoticeable in obstretncs since the fetus is often quite active and 
good images for measurements are rrot easy to obtain m any 
defined order The analysis packages are also available during 
video tape playback so free-form inpu1 helps to etfrninate un- 
necessary tape searching that would be required to find the 
current measurement asked for by the walk-through method 

The operator gains access to the quantitative measurement 
packages by pressing tine ANALYSIS key on the 77020A's 
keyboard. Upon entering ANALYSIS mode, a menu is presented 
to the operator which associates a measurement number with 
the measurement name and identifies the special function keys 
(see Fig. 1 ). An operator prompt message, selB, is also printed 
on the screen To store a measurement as analysis package 
data, the operator positions the system calipers by pressing the 
CAUP 2 key and using the joystick, and then typing the number 
associated with the measurement. When the numerical key is 
typed, the system stores the value indicated by the calipers into 
memory as one value for that parameter and the value is dis- 
played on the screen below the caliper values, Several values 
for each measurement can be stored and the average of these 
is used by the system when performing the analysis package 
oalculations. 

After measurements are made, the operator presses the 
APPROVE key. The display is cleared and the operator is pre- 
sented with the analysis package edit report (see Fig 2) This 
report contains patient identification, time and date, patient his- 
tory, all stored measurements and their averages, the configured 
calculations, and comments. A menu of commands gives the 
operator the ability to modify the report as needed. 

Two sections of the edit report are reserved for qualitative 
information about the patient and the examination. The patient 
history section of the report contains configurable prompt head- 
ings followed by enough space to insert the appropriate informa- 
tion. This section contains relevant patient information such as 
age, height, and weight. The edit history menu command is 
used to add or modify the values associated with each heading. 
The COMMENTS menu command can be used to enter three Sines 
of text summarizing the findings for this particular patient. 

The quantitative information is summarized in the measurement 
and calculation sections of the edit report. These sections are 



ment with the joystick (see Fig, 1 on page 8}, The paddle 
switch changes the size of a displayed ellipse and allows 
measurement of an area or volume (the ellipse rotated 
around the cursor axis). 

The CAUPer 2 key freezes the current set of measurements 
and enables a second set. After the second set is enabled, 
the CAUP 2 button toggles between the two sets, leaving 
both lines or ellipses on the screen, ANALYSIS allows image 
analysis using the internal obstetric or cardiac analysis 
software (see article on page 8). SET TIME allows the internal 
calendar clock to be set to the current date and time. This 
data, along with all other setup data, is stored during 
power-down conditions in nonvolatile RAM, 

SECTOR SIZE switches the size of the sector on the monitor 
from normal to large, allowing the diagnostician to see 
more detail in the image, CONF puts the system into a con- 



dependent upon one another since the calculated values are 
computed biy the ultrasound system from the appropriate av^* 
age measufement values. Wtien the operator selects the EDU 
MEASURE command, the operator can add. rr^xjrfy, or delete 
measurements as needed- Any ediied measurements are 
marked with an asterisk to indicate tha! the value was manually 
entered and not measured by trie system When the operator <s 
finished editing measurements, the measurement averages and 
calculations are recomputed and redisplayed 

The operator may wish lo edit the calculation values tor several 
reasons. CalculatJons that are rrretevant for the cunent patient 
can be deleted The operator may wish to use a calculation 
fomiula that is unavailable in the obstetric package; editing allows 
the operator to compute the value and type it in ihe appropriate 
space manually In the cardiac pacl<:age, a different volume ah 
gofithm can be selected Jor the calculation if the operator feels 
that it is more appropriate for the current patient than the config- 
ured one. Selecting the edjt calculations menu command 
allows the operator to edit the calculation values or change a 
volume method. Any edited values are marked to indicate that 
they were edited and not computed by the system. 

After any necessary editing is completed, the operator has 
tv^/o choices. If it is fett that more measurements should be made, 
the operator can return to reaMime imaging or tape playback by 
choosing the image menu command. All stored measurements. 
history, and comments are retained and any new measurements 
are appended to the report if the operator is satisfied with the 
contents of the edit report, the operator selects the final report 
commartd from the edit report menu The screen is updated with 
the final version of the reF>0Jl. The final report differs from the 
edit report in that only the average values tor each of the measure- 
ments are displayed and it is reformatted slightly for aesthetic 
reasons. No edits can be made lo this report. Any of the system 
hard -copy devices can be used to make copies of the final report 
for filing purposes. 

The analysis packages are compfemented by the system cOn- 
FiGURATtON command, This command ailows the operator to 
choose certain system parameters and anaiysis package fea- 
tures Menus are used to configure the anaiysis package fea- 
tures. One menu allows the operator to select from a list of avail- 
able calculations the ones the system will perform aulomatically- 
Another menu allows the operator to configure the patient history 
headings as appropriate for the particular department using the 
analysis package The configuration data for each analysis pack- 
age is stored separately in system memory to ailow for differences 
between the two packages 



figuration mode, which allows the doctor to customize the 
system for the hospital's needs. 

System Description 

Fig. 2 shows a block diagram of the 77020A, Based upon 
the particular transducer, the depth setting, and whether 
M-mode is selected, the puiser and timing section of the 
scanner generates 64 different transmit pulses » each at the 
proper frequency and phase to steer the acoustic beam in 
the desired direction. These are amplified by the drivers 
to the level determined by the transmit control, and trans- 
duced into acQUStk energy by the transducer array. 

The echoes detected by the transducer have about 115 
dB of dynamic range. Displaying this information to the 
observer who. at best, has a limited visual dynamic range 
of 25 dB, requires much processing. First, the TGC settings 
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Fig, 2. Block diagram of the 77020 A Ultrasound imaging System. 



applied to the 64 preamplifiers window an appropriate 
55-dB segment of the retiirniiig signals. This data is then 
combined in the phasing delays to generate the video sig- 
naL Based upon the operator settings of the COMPRESS, SW 
GAIN, and PRE-PROC controls, this video signal is processed 
and sent to the stripchart recording interface circuits in 
the controller for paper-trace hard copy and to the display 
subsystem. 

Before it is presented to the display subsystem, this data 
can be intercepted by an HP 1000 Computer for analysis 
or tissue identification studies. Allowing this to happen 
before scan conversion of the acoustic data gives the system 
the capabilities of real-time data capture without any of 
the data reduction problems that conversion to standard 
NTSC/CCIR video generates. The raw RF acoustic data can 
also be captured. 

This data now enters the digital scan converter where it 
is converted from the acoustic data's polar coordinates to 
raster-output Cartesian coordinates using a special conver- 
sion algorithm to remap the data (see article on page 30]. 
Stored in a 480x640 five-bit display memory and 
postprocessed (gray-scale mapped] as selected by the 
operator, the data is summed with graphic/character data 
from the bit map and converted to video for display on the 
high-speed display in the controller, storage on the VCR 
for later playback, transmission to any standard video aux- 
iliary device, or transmission to the video camera for haid- 
copy output. 



The high-speed display in the 7 7900 A Controller was 
especially developed for 60-Hz noninterlaced (flicker-free) 
video, gamma -corrected for the human eye. The controller 
is the node for both control and peripheral data conversion. 
All controls hard copy, and storage data is routed through 
this instrument to allow simple, cost-effective updates. The 
control panels are also pari of the controller. 

To allow full gray -scale control of the video playback 
from the VCR. the image is converted from analog to digital 
data and read back into the memory. Audio data such as 
the doctor's verbal comments or diagnostic heartsound data 
is also stored and played back through this system. 

Stripchart recording of gray- scale data requires a fiberop- 
tic oscillographic recorder. This normally is used for M- 
mode hard-copy output. When the operator wishes to add 
a full-frame image to the M -mode stripchart recording, data 
can be routed via the internal HP -IB (IEEE 488) from the 
digital scan converter memory through the controller to 
the stripchart recorder. This data is also gamma-corrected 
for the dry-silver media used by the recorder so that the 
hard-copy image matches the display. The operator also 
has the option of adding physiological data such as elec- 
trocardiograms to the image or hard copy for timing and 
additional diagnostic information. 
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A Physician's View of Echocardiographic 
Imaging 



by Richard L. Popp, M.D. 



A FEW PHYSICIANS ARE INTERESTED primarily in 
ultrasound, noninvasive imaging, or other applica- 
tions of technology to medicme. But most physicians 
and health care professionals want medical infonmation 
that will help them care for patients* This majority becomes 
interested in an imaging device if it provides new data not 
availahle by other means, more complete or more accurate 
data than competing methods, advantages such as obtain- 
ing data at less risk or less cost to the patient, or more con- 
venience to the patient and/or physician. Ultrasound imag- 
ing h attractive because it can sat i sfy all these conditions. 

One contribution of ultrasound imaging is the ability to 
*'see'* the heart or other structures inside the body without 
doing anything potentially dangerous to the patient. 
Another feature is the detail with which the anatomy of 
soft tissues, such as heart valves and heart walls, are dis- 
played. Heart valves and walls move rapidly and no other 
widely available imaging method can show both structure 
and motion with such detail as two-dimensional ul- 
trasound imaging (2-D echocardiography]. Strictly speak- 
ing, I can't say that new conditions or diseases have been 
discovered by echocardiography, but this method has given 
us new information to help explain things we have only 
observed about patients in the past. For example, some 
patients develop heart failure after heart valve replacement 
surgery. Echocardiography is easily performed in the inten- 
sive care unit as well as in the laboratory and allows us to 
decide if the heart muscle or valve is malfunctioning. We 
can also tell if the heart is compressed by fluid accumulat- 
ing within the chest or if the valve itself has pulled loose. 
In this example echocardiography is done rapidly and pro- 
vides accurate information conveniently, 

Cllnicaf Use of Ultrasound Imaging 

In general, echocardiogrtiphy has made an important con- 
tribution to patient care because it is noninvasive [nothing 
need be injected into the body), it can be done without 
hospitalizing the patient [and so the cost is kept down), 
there is no ionizing radiation involved, and there is no 
known risk to the patient from performing the procedure. 
It offers improved differentiation of soft tissues which look 
the same in X-ray images and because of all of the above, 
it is quite repeatable. The repeatability translates to the 
ability to follow patients and learn the natural history of 
their condition with and without interventions to change 
that natural history. Jf the physician simply wants to make 
a diagnosis of a suspected condition, the physician might 
often use echocardiography now» whereas a cardiac 
catheterization, with attendant risk and cost, would have 
been the method of choice in the past. Certainly the 
echocardiogram is superior to any other technique for the 
diagnosis of small masses within the heart such as cardiac 



tumors and bacterial growths on the heart valves, and can 
better display fluid accumulations around the heart. 

Making the statement that echocardiography is superior 
to other techniques has required a long process of testing. 
From a historical perspective one must realize that pulsed 
echo techniques date from the mid 1950s and high- volume 
use of echocardiography in the United States, and sub- 
sequently in the world, has evolved since approximately 
1967, During the early part of this period, M-mode echocar- 
diography was the only form available. This mode provides 
information on the dynamics of cardiac motion from only 
a single area of the heart at any given time since the ul- 
trasound beam must be reoriented to encounter another 
part of the heart for study. The need to expand the display 
to provide an overall picture of the heart in anatomic format 
was appreciated, and 2-D echocardiography was invented 
in the early 1970s- This was quite rudimentary until the 
mid 1970s and the modem era of 2-D echocardiography 
has existed only the past five years or so. Although much 
of the data that is useful clinically is available through 
M-mode echocardiography for those who have learned how 
to interpret it, 2-D echocardiography makes the same data 
more comprehensible to more physicians and provides new 
data not available by M-mode. 

Comparison of M-Mode and Two-Dimenslonal Methods 

The extension of M-mode echocardiography to 2-D ul- 
trasound imaging usually means increased cost. For this 
reason we must explore the relative value of the old and 
new methods. The high sampling frequency of M-mode 
(about 1 kHz] compared with 2-D echocardiographic im- 
ages [about 30 Hz) means that very fast cardiac movements 
are tracked better with M-mode. For example, high-fre- 
quency vibrations of valve structures indicating abnormal 
flow are usually more obvious in M-mode than in 2-D 
images, and in fact may be imperceptible in 2-D studies. 
The subtleties of motion of individual echoes may be more 
easily analyzed by M-mode because 2-D images also seem 
to present information "overload" to the observ^er. That is, 
so much data is presented on a 2-D image that the observer 
may have trouble concentrating on a single echo among all 
of the signals making up the 2-D image. Very slight separa- 
tions of the pericardial sac surrounding the heart seem 
more clearly defined, for example, using the limited sam- 
pling area of the M-mode technique. The human eye/brain 
combination has difficulty keeping track of both cardiac 
motion and the time base provided by the electrocardio- 
grara* Thos> it may be easier to match events on the M-mode 
records when compared with 2-D images. 

On the other handt the expanded anatomic information 
presented by 2-D echocardiography has revolutionized the 
diagnosis of congenital heart disease, helping us to under- 
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stand the relative size of the heart chambers and the number 
of heart valves that are properly and improperly formed. 
The alteration in cardmc wall motion relative to more nor- 
mal motion among segments of the heart is used to recog- 
nize and quantify the effects of coronary arter^^ disease on 
the heart. This is done almost exclusively with 2-D echocar- 
diography. The abnormal bulge (aneurysm) that is the result 
of heart attack In many patients is recognizable by 2-D 
methods and is not perceptible hy most M-mode methods. 

In fact, the combination of 2-D and M-mode methods 
often is optimal for understanding each type of record and 
best serves the patients Because virtually all 2-D instru- 
ments are also capable of providing M-mode recordings, 
we need not pit M-mode against 2-D methods for diagnosis. 
It is more appropriate to compare the limited ability of 
M-mode methods alone versus the combination of 2-D im- 
ages and M-mode recording for the comprehensive analysis 
of each case. 

Although I have mentioned that current echocardlo- 
graphic methods are extremely useful and that they are a 
relatively recent development, one may properly ask why 
physicians are not rushing out even more quickly to equip 
themselves with 2-D echocardiographic instruments? To 
answer this question we should first become familiar with 
the details of the information sought from the images 
created* With respect to the heart, the images may show 
abnormal structure such as a congenital malfunction, or 
distortion of the originally normal anatomy by disease pro- 
cesses such as rheumatic heart disease. The presence of 
structures that are not normally present, such as blood 
clots, tumors, or masses of bacteria within the heart, are 
seen and recognized as abnormal features (Fig. 1). The "bet- 
ter" the quality of the image produced by the instrument , 
the better our ability to recognize subtle abnormalities of 
structure, distortion of anatomy, or abnormal new struc- 



tures. Anatomy also is important when one looks at the 
secondary effects of disease such as enlargement of heart 
chambers and altered shape of these chambers. 

In addition to simply looking at the size and shape of 
the heart, the dynamic parameters of cardiac motion are 
of great interest to the cardiologist. Some abnormal patterns 
of motion alone represent disease processes such as those 
associated with rheumatic mitral valve disease. The 
amount the heart w^alls move tells us about the volume of 
blood moving through the heart and the rate of motion 
reflects the rate of volume change or flow. In addition to 
measuring the size and shape of the organ and observing 
the dynamics of the moving cardiac structures, one also 
can notice normal or abnormal appearance of the tissue 
being displayed. These alternations in tissue appearance, 
or tissue texture, are of special interest to the radiologist 
imaging abdominal organs, for example, but also are of 
great intj^rest to the cardiologist. Currently we are able to 
recognizee calcification of tissues such as the diseased aortic 
or mitral valve, scars that occur within the heart after heart 
attack, and infiltration of the heart with abnormal proteins 
such as amyloid. The ability to see all three types of infor- 
mation (anatomy, function, tissue integrity) is w^hat deter- 
mines a ''good" system. 

Equipment Considerations 

Observing the size and shape of the heart, or other organs, 
means looking at measurements that depend on the resolu- 
tion of the instrument. As resolution improves and the 
image becomes more clearly defined, it becomes easier for 
the user to outline the heart to measure its size and extrapo- 
late this size to volume. A clear outline throughout the 
cardiac cycle leads to accuracy in measuring dynamic pa- 
rameters. Differentiating very low-intensity signals from 
the background noise on the image is also important be- 




Ffg* 1. An echocardiographic 

image of the heart taken with the 
trar^sducer on the chest waH in the 
parasternat position and ofiented 
parailei to the short axis of the ieft 
ventricle (Iv). The arrow indicates 
an abnormal mass within the right 
ventricle {r/} attached to the inter- 
ventrjcuiar septum. The mass has 
visual features that make tt appear 
different from the myocardium 
(heart muscle). The mass is a 
malignant tumor from a distant site 
which has traveled to and grown 
in the heart. Note the echo- free 
(black) ventricular chambers 
compared to the low amplitude 
echoes within the myocardium 
and the tumor. Wail thickness can 
be judged from this image even 
in the absence of specular echoes 
from the blood-muscle interface^ 
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cause it is the combination of very low-intensity and ver\' 
high -in tens it>^ signals that produces the texture one uses 
to assess the quaiit>^ of the tissue (Fig. 1), 

As the various aspects of these considerations become 
optimized, ultrasound images become more and moresimi* 
lar to optical images of the organ being studied. Very clear 
pictures of the structures and lack of spurious instrument- 
generated echoes please the user (Fig. 2). Other features 
such as gray-scale mapping and processing schemes can 
alter the aesthetics of the picture, but these merely satisf>^ 
individual user preference, by and large. Optimizing the 
system also includes maidng it easy to use. This includes 
not only putting the minimum number of controls on the 
instruments but also having these controls in the best pos- 
sible location for use and making them easy to activate. 
Other considerations come from clinical use on hos- 
pitalized patients and these include the size of the trans- 
ducer and the portability of the instrument. 

More resolution is obtained by using a maximum aper- 
ture transducer, hut the space between the ribs must be 
used for the sound entering the body and this means that 
a transducer face larger than 20 to 25 mm gives difficult>^ 
vv^hen applying it to the chest wall. The nee^d to manipulate 
the transducer on the chest wall, tilting it and rotating it, 
produces special consi derations for those designing the 
transducer case. 

Finally, the acceptance of echocardiography has given 
rise to a need to use it to make clinical decisions on very 
ill patients, sometimes in life and death situations. As a 
result one often needs to transport the instrument to the 
intensive care unit, squeeze it between the patient^s bed 
and the life support equipment, and rapidly record infor- 
mation in the midst of feverish activity around the patient. 
This requires a small, light, maneuverable, and easily used 
system with self-contained recording equipment. 

The historical development of mechanical 2-D scanners 
and electronic phased-array 2-D scanners has led to some 
generalities about the equipment that is prevalent in the 
user public. The original advantage, the wide angle of view, 
of phased-array instruments over mechanical scanners is 
no longer present. The ability to record simultaneously 
M-mode data and 2-D images is mainly the property of 
phased-array scanners, although some mechanical instru- 
ments now have one transducer for 2-D imaging and a 
second transducer for simultaneous M-mode recording. 
Mechanical scanners are usually less expensive than 
phased-array scanners, although the gap is variable and 
narrowing. Mechanical scanners have a reputation for a 
less cluttered, low-background -noise image compared to 
phased-array devices. The phased-array systems originally 
had significant problems with spurious echoes contributed 
by grating lobes and generally had more complex controls. 
Phased-array instruments of the generation characterized 
by the HP 77020 A generally are more sensitive when judged 
in the clinical sphere and compared with mechanical scan- 
ners. That is, the ability to get usable information on a very 
high percentage of patients seems better with these instru- 
ments than with most mechanical scanners. In general, the 
clean, clear images and "beautifur' pictures produced by 
the best mechanical scanners in patients that are easy to 
image have had a slight edge over the phased-array images 



for these patients. However, with even slight increases in 
imaging difficulty caused by biologicaJ and technical fac- 
tors, the phased-array devices usuaOy have been able to 
acquire clinically useful data when the mechanical scan- 
ners have not been up to the job. 

General Thoughts on Echocardiogra|ihic Testing 

To place the current application of echocardiography in 
proper perspective for its evaluation in the cHnical setting, 
it is worthwhile to review the methods by which such 
appUcations are developed. With the invention of echocar- 
diography, as well as the subsequent improvements cul- 
minating in the most advanced equipment, a general pro- 
cess applies. First, patients with well-confirmed or known 
disease are selected for testing with the new equipment to 
see the result. In this group of patients with 100% preva- 
lence of the condition, the patterns of abnormality are de- 
fined and can be learned as examples of the ''classic case," 
Assuming that there are characteristic abnormalities on the 
test accurately indicating a certain condition for the pa- 
tients known to have the condition, then the method is 
applied to patients who are suspected of having the disease 
or process under discussion. In these patients the preva- 
lence of the condition may be very low or ver^' high, gen- 
erally depending on the degree of rarity of the condition 
and also depending on the clinical acumen of the physician 
who has taken the history, done the physical examination, 
etc, J before considering the new test. 

For physicians to use echocardiography properly in vari- 
ous clinical settings, it is important to have data on the 
sensitivity and specificity of the test for the condition under 
discussion, to know the prevalence of the condition in the 
population, and to judge the probability of this individual 
patient having the condition sought. Learning to rely on 
the test when it is positive and/or negative requires a sig- 
nificant data base obtained from doing the test in patients 
with a well-defined presence or absence of the condition 
as validated by other accepted techniques. Such informa- 
tion is needed before we are able to say either that a con- 
dition truly is present when the test is positive, or that we 
can exclude the condition when the test is negative. If a 
physician is trying to confirm a strong clinical suspicion 
that the disease is present, then the physician wants a test 
that confirms the presence of disease when the test is posi- 
tive or abnormal. Echocardiography does this extremely 
well. Conversely, if the clinician is trying to reduce the 
possible diagnoses in a given patient by progressively rul- 
ing out specific disease states, this process requires the 
clinician to exclude the disease confidently. In general one 
cannot exclude disease processes well with echocardio- 
graphic images. For example, when a large mass is present 
within the heart, we see it well. When we are trying to 
exclude the possibility of a mass, and none is demonstrated 
on the image, we are still left with the question whether 
there truly is no mass or whether it is so small that it is 
beyond the resolving capability of the instrument. 

A related issue develops from application of echocar- 
diography in some conditions where new information is 
uniquely available. For example, the presence of small 
masses, such as groups of bacteria on a heart valve, may 
be recognized with echocardiography when no other cur- 
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Fig, 2. r/vee sa/'/a/ sections (images) parallel to the long axis of the left ventricle (iV) from the 
patient shown tn Fig. 1 The mass within the right ventricle and right venthcufar outflow tract 
(RVOT) contains very small blood vessels (arrovvs) that can be resolved with this imaging 
system. Sena/ sections of this type allow the echocardiographer to create a mental image of 
the spatial relations of the heart. (A0= aorta.) 



rently available technique can do this. If this new test is 
superior to the old test, then the new test seems valuable 
intuitively, but it is very difficult to prove this from a sci- 
entific standpoint. If bacteria are present in the patient's 
bloodstream and the echocardiogram does not show a mass 
within the heart, it is impossible to say whether this patient 
has an infection on the heart valve that is not shown by 
the echocardiogram, or whether in fact the bacteria in the 
bloodstream come from some other source within the body. 
This kind of consideration makes the establishment of sen- 
sitivity and specificity values for a new test very difficult. 

We may consider using a relatively attractive noninva- 
sive test in the screening of asymptomatic populations to 
look for disease processes. Again, depending on the preva- 
lence within the population of the abnormalities sought, 
we can determine if using this test is worthwhile only 
when we know the sensitivity and specificity of the test 
for conditions we seek. Echocardiography possibly could 
be used for screening for some conditions, but at present 
the data base is not sufficient to help us determine whether 
this is worthwhile. We are in a phase of evolution in the 
technique and many of the expanded uses of the technique 
await data of this type. We can determine whether echocar- 
diography is worth the cost only if we know both the quan- 
titative accuracy of the method and when we know the 
impact of such information on long-term patient manage- 
ment. An example comes from using ultrasound imaging 
to help manage patients after heart attack. Most clinicians 
think they would like to know how much heart muscle 
remains, how well it contracts, and whether an aneurysm 
or abnormal scar is present after a heart attack. To get such 
information now, an angiogram or radionuclide study must 
be performed. If the patient can have echocardiography 
without radiation hazard and with reduced cost compared 
to these methods, and if the test is accurate, we can then 
design studies to learn if this information affects patient 
management and prognosis. 



Training 

The studies showing that echocardiography is very use- 
ful clinically and very accurate quantitatively have been 
done predominately by those with great skill and experi- 
ence in the technique. Primary training in performance of 
echocardiography and interpretation of the results has been 
included in formal cardiology, radiology, and ultrasound 
training programs only within the past ten years* Recogni- 
tion of the value of the method has created a tremendous 
need for training. Those people currently responsible for 
training people in echocardiography generally agree that a 
significant period of time is necessary to perfect the skills 
needed to perform the technique properly for clinical deci- 
sion making. Guidelines have been established by the 
American Society of Echocardiography, and endorsed by 
other organizations, for both physician and technician 
training in both M-mode and 2-D echocardiography. 

Between three and six months of physician training and 
experience is needed for independent interpretation of M- 
mode and 2-D echocardiography. This time commitment 
makes training quite difficult for the practitioner who may 
not be able to suspend a medical practice. Several proposals 
are now under consideration for dealing with this problem, 
but lack of proper training opportunities means echocar- 
diography may be performed suboptimally in some settings 
and patients who would benefit from its optimal use may 
be denied this clinically important method. 
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An Acoustic Transducer Array for Medical 
Imaging — Part I 



by John D. Larson III 



THE ACOUSTIC TRANSDUCER ARRAY is the analog 
front end of an ultrasound imaging system. It pro- 
\ades a large number of independent channels, trans- 
duces electric signals to acousdc pressure, and generates 
sufficient aconstic energy' to Illuminate the various struc- 
tures in the human body. In turn, it converts the weak 
returning acoustic echoes to a set of electric signals which 
can be processed into an image. The key transducer require- 
ments are to: 

■ Generate 10 to 100 mW/cm^ of acoustic power from 
reasonable input voltages 

■ Provide good signal- to-noise ratio 

■ Have low cross- coup ling between elements 

« Generate short acoustic pulses, =^2 ^s in duration 

■ Couple acoustic energy to the patient efficiently 

■ Damp the aconstic backwave 

■ Achieve broad angular coverage 

■ Suppress un desired vibration modes 
m Be lightweight and handheld- 

Although a variety of techniques for generating and de- 
tecting sound are available^ piezoelectric transducers are 
preferable because they can both generate and receive 
acoustic waves » and they are several orders of magnitude 
more sensitive than other possible transducers.^ 

To construct an array from typical ceramic piezoelectric 
materials, a number of fabrication and performance prob- 
lems must be solved. This article discusses the fabrication 
of a transducer array for medical imaging and some of the 
basic principles and design constraints that must be consi- 
dered. The following article on page 22 describes the appli- 
cation of these considerations to the design of the transduc- 
ers for the HP 77020A Ultrasound Imaging System. 



Couptmg to Patient and Short Pulse Generation 

A pulsed 2. 5 -MHz piezoelectric transducer with no 
matching or damping layers rings and produces a long 
pulse on the order of 10 /iS or more. Because range resolu* 
tion is inversely proportional to pulse lengthy this reso- 
nance is unacceptable. Placing the transducer in contact 
with the patient neither damps this resonance greatly nor 
couples energi^ very efficiently because there is more than 
a 20:1 acoustic mismatch between currently available 
piezoelectric materials and the human body. 

To obtain a shorter pulse that is more effectively coupled 
to the patient, several techniques can be used singly or in 
combination. The goal is to reduce the 20-to-30'dB inser- 
tion loss caused by the large impedance mismatch and to 
produce pulses short enough to allow 1-mm range resolu- 
tion- 

Ceramic Resonator in Direct Patient Contact The loaded 
quality factor Ql of a ceramic resonator in contact with a 
mismatched impedance [Fig, la) and possessing its own 
internal mechanical Q, Q^, is given by:^ 



Ql= -^ 



1 + r^A 
t-r^A 



where A = exp( - 2tt/Q^] and r is the reflection coefficient. 

For a typical lead-zirconate-titanate (PZT) ceramic, 
r = 0.8S> Q^-70, and Q^^ is 18-4. The time response is a 
damped sinusoid of the form 

exp(-TTf^,t/Qj^) stn[2irf^t) 
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Fig. 1. Schemes to obtain short 
pulses of uHrasound from acoustic 
transducers (a) Undamped pi- 
ezoelectnc transducer placed di- 
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coupling between the patient and 
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(a) using transwrsai fifter pro- 
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For a resonator with fo = 2,5 MHz^ it takes 8.1 fxs for the 
response to drop 30 dB, and the resulting range resolution 
{6 mm) is unacceptable. 

Matched, Lossy Backing. If an acoustically matched and 
absorbent backing is added to the rear of the transducer 
(Fig, lb), the effective Qi, is lowered to ^2 and the pulse 
length is reduced to ^1 ^s. The acoustic energy propagated 
into the backing, or backwave, is lost by acoustic attenua- 
tion and no secondary resonances are set up. For example, 
for lead metaniobate ceramic material operating at 2,5 MHz 
and transmitting ultrasound into human tissue* the pulse 
duration is L3 /as. The insertion loss, including transduc- 
tion efficiency and the impedance mismatch of ceramic to 
tissue is -28 dB, 

This technique clearly improves the puise response, al- 
though most of the power (>99%) is lost in the backing. It 
also has the advantage of minimizing the number of acous- 
tic bonds and acoustic parts. 

Matching Layers Between Transducer and Patient. In this 
realization [Fig. Ic), one or more matching layers are used 
to improve the transmission of acoustic energy into the 
body. The effective Q^ of the resonator is reduced, so a 
shorter pulse of ultrasound can be obtained. For example, 
a typical PZT ceramic with two matching layers of appro- 
priate acoustic impedances interposed between the 
ceramic and the patient has a -4-dB insertion loss and a 
2.1-^s pulse duration. This represents a large improvement 
in transduction efficiency with some loss of pulse length 
compared to the case shown in Fig, lb. 
Inverse Filter. As mentioned earlier, the pulse excitation 
of an unbacked transducer (Fig, id] yields a damped 
sinusoidal waveform, A compression filter might be built 
that takes a sinusoidal waveform and yields a pulse output. 
One such realization is an inverse filter.'' In principle, sub- 
stantial pulse compression can be obtained. In practice, 
the finite noise in the received signal limits the filtering 
process and the compression is limited to reduction factors 
of 2 to 3 in pulse length."* 

The techniques illustrated in Fig. lb and Fig. Ic are used 
in the 7702GA Ultrasound Imaging System because of good 
pulse response and decreased insertion loss. These attri- 
botes lead to lower drive voltages, higher reception sen- 
sitivity, and pulse lengths commensurate with the desired 
range resolution. 

Packaging Requirements and Constraints 

A basic requirement of a phased-array transducer is that 
the elements be placed on half -wavelength [X/2) centers to 
sample the aperture properly. Here \ is measured in the 
subject material For the 1.540-2.5 -cm apertures used for 
medical imaging, 60 to 80 elements are required. In turn, 
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Fig, 3, Basic acoustic transduoef array configuration. 

the ceramic transducer must be a ha If- wavelength (A/2) 
thick to operate properly. Here A is measured in the 
ceramic. Thus, the elements have a width-to-thickness ratio 
directly proportional to the velocity of sound in the 
piezoelectric ceramic and inversely proportional to the ve- 
locity in tissue. This ratio is about 3:1, resulting in tall 
narrow elements. 

Fig. 2 illustrates the possible lead connection schemes* 
An independent lead to each element is required. In Fig, 
2a> the high-voltage electrodes (=^200V pulse present) are 
toward the patient, causing a safety hazard. To overcome 
this, the high -voltage electrodes can be located on the back- 
ing side of the element if the backing material is an insulator 
(Fig, 2b). 

The preferred structure is illustrated in Fig. 2c. Here an 
insulator is placed between the backing and the ceramic. 
A metal foil is bonded over the elements to form the ground 
connection, protect the patient, and keep grease out of the 
kerfs between elements (the elements are cut apart to sup- 
press cross-coupling). 

An important modification is to use a second layer of 
piezoelectric ceramic as the insulating layer. This gives 
suppression of the unwanted mass-spring mode and in- 
creases the acoustic energj^ radiated by the transducer. 

Array Element Design 

The acoustic arrays used in the 77020A's handheld trans- 
ducers are composed of a large number of identical ele- 
ments. These elements must be protected from the outside 
world of grease and dirt, and must be packaged with a 
suitable lens and case to provide a useful probe. 
Driving Point Considerations. Each array element appears 
electrically to be a frequency- dependent, lossy capacitor 
with capacitance Cj^^^tjWL/H where e^ is the clamped 
dielectric constant and W, H, and L are defined in Fig. 3. 
As will be seen later, these dimensions are constrained by 
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the need to control spurious modes and grating lobes, and 
to obtain a narrow beam in tlie elevation direction. 

A compromise between high-frequency, shallow-pene- 
tration operation on one hand and low-frequency, deep- 
penetration operation on the other leads to choosing a 
center frequency of 2.5 MHz, a frequency which is assumed 
in the examples to follow. 

To S3tisf>' the above conditions, the dimensions L = 250 
^m. period P = 320 imu, W = io mm. and H = 600 fxm for 
PZT-5H ceramic material were chosen, and Co =55 pF. 
Transmit Considerations. Exposure of mammalian tissue 
to average power densities of 100 mW/cm" or less has been 
found to cause no significant biological effecls,^'^ The de- 
sign goal has been to stay below this power leveL Detailed 
calculations show that for PZT-5H, a single-layer array ele- 
ment can achieve 100 mW/cm^ for a 200:1 duty cycle with 
216 volts peak'to-pealt drive, A double-layer element devel- 
ops the same acoustic power output at 130 volts peak*to- 
peak. The use of acoustic matching layers can allow more 
efficient power generation. For example, a single matching 
layer used in conjunction with a single PZT layer could 
result in a transduction efficiency increase of 8 dB, or 
around 90 volts peak-to-peak to achieve 100 mW/cm^. 

Another piezoelectric material of interest is lead 
metaniobate* It is similar to PZT-5H, but its dielectric con* 
stant is only one-sixth as great and it has a lower elec- 
troacoustic coupling coefficient. Its chief attraction is a low 
spurious mode level. For a 200 ;1 duty cycle and 100 mW/ 
cm^ power level, a single lead metaniobate layer requires 
540 volts peak-to- peak drive while a double layer requires 
360 volts peak-to-peak. The high drive voltage makes it 
less suitable for use with solid-state circuitry. 

Since the acoustic arrays are operated in a pulsed mode 
to achieve good range resolution, wideband operation be- 
comes a prime consideration. Factors influencing this in- 
clude frequency-sensitive driving-point impedance, epoxy 
bond thickness, and frequency-dependent transduction ef- 
ficiency. Other vibration modes can cause narrowband. 
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extended -time- duration response. 

Using a circuit model for an acoustic transducer, the 
driving-point impedance of a single-layer PZT array ele- 
ment can be calculated as given in Fig. 4. The magnitude 
decreases with increasing frequency while the real part 
reaches a peak around the 2.5'MHz center frequency. Fur- 
thermore, by comparing measured data to the theoretical 
data, it is apparent that there are high-Q modes of vibration 
present which the simple circuit model does not predict. 
Epoxy Bond Effects. The transducer response is quite de- 
pendent on the quality of the acoustic properties of the 
bonded array. Chief among these are the thicknesses of the 
epoxy-bond layers used to cement the PZT to the backing, 
the surface layers present on the backing, and the degree 
to which the backing impedance is the same as the PZT 
acoustic impedance. 

The epoxy bond has an acoustic impedance of 3.2X10^ 
ray Is. compared to 25x10*^ ray Is for PZT-5H. This 8:1 mis- 
match dictates that the bond be very thin to prevent narrow- 
ing of the passband. Fig, 5 shows how the frequency re- 
sponse of the transducer is changed as the bonds vary from 
1 pLin thick to 8 /im thick. In turn, the bandwidth varies 
from 74% to 31% as this occurs. To ensure good response, 
bond thicknesses less than 2 ^m are required. 
Acoustic Modes of Narrow Elements. In an acoustic imag- 
ing array the elements should have a single mode of vibra- 
tioOf namely a piston-like motion of the front surface. By 
making the elements sufficiently narrow, and reducing 
cross-couplingj the acoustic intensity should be uniform 
at various angles from the normal (see Hg. 3), that is, om- 
nidirectional. For illumination out to angles where the re* 
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sponse is no more than 3 dB down from thai at normal 
incidence, the element width L should be ^OA4KIsin0^^Q 
which for ±45° coverage and 2.5-MHz operation, yields 

Fabrication constraints force the actual elements to be 
even narrower With the element thickness H comparable 
to Lniax* additional modes can be expected to exist at fre- 
quencies near the desired thickness mode frequency. Such 
vibrations include the dilatational or breathing mode, the 
mass-spring oscillation mode, Lamb wave on a foil, and 
the Rayleigh surface- wave mode. 

Onoe and Tiersten^ give a good account of the thickness 
and dilatation modes. Defining a configuration ratio G — L/ 
H, it can be shown that these modes are close in frequency 
and strongly coupled when G-L For G<<1, the modes 
become widely separated in frequency. 

Fig. 6 illustrates these two modes for 760-/im-thick PZT- 
5H at G = 0.1 to 1.3. The thickness mode frequency is 
essentially constant at 2 MHz while the dilatational mode 
frequency is inversely proportional to width L. 

Another major acoustic vibration mode is termed the 
low-frequency mode^ because of its manifestation at fre- 
quencies around 0.8 to 1 MHz, as compared to the main 
thickness mode at 2 MHz. It occurs as a fairly high-Q, 
slowly damped resonance which requires 40 to 50 ju.s to 
decay below the system noise level. The origin of the mode 
appears to be a harmonic oscillation of the element moving 
against the spring formed by tlie backing material. The 
frequency of this mass -spring mode can be given by 
^m-a - VeffVl/LH, where v^fj = vc^pp^, c^ is the elastic con- 
stant of the backing, and ppzT 1^ the density of the piezoelec- 
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trie element. 

The main effect of the mass -spring mode on transducer 
performance in an acoustic imaging situation is to add 
clutter to the image for the range from to 7 cm because 
the mode is highly excited by the transmit drive pulse. 

A thin metal foil is bonded to the top surface of the array 
to seal the element interstices. This foil can support a 
guided wave known as the Lamb wave,** Two subgroups 
can be identified depending on particle motion of the foil. 
The symmetric modes show particle motion where both 
sides of the foil move away from each other or towards 
each other* They are relatively nondispersive, that is^ the 
acoustic velocity does not var\' greatly with frequency. The 
asymmetric modes exhibit flexural propagation with ve- 
locities approaching zero as the frequency or plate thick- 
ness goes to zero, hi metal foils, the symmetric mode has 
a velocity of about 4x10^ cm/s, while the asymmetric mode 
varies from to 2 x 10^ cm/s. The excitation of these modes 
comes from the motion of a driving element, and is man- 
ifested mainly by cross-coupling between elements. 

The last major mode that exists on the array is a Rayleigh 
surface wave.^° This propagates on the backing with vol- 
ocities in the range of 1 to 2 x 10^ cm/s. The particle motion 
is a retrograde ellipse composed of components along the 
surface and perpendicular to it. This motion decays expo- 
nentially with depth from the surface. This mode is driven 
by the motion of an element and is another cause of cross- 
coupling between elements. 

Mode Suppression. Of the modes discussed, the thickness 
mode is desired and all others are undesired. The dilata- 
tional mode is not suppressed, but by choosing narrow 
elements, the mode is moved to a higher frequency where 
it is conveniently filtered out. 

Of the remaining modes, the mass-spring mode is most 
troublesome. It is strongly excited and has high Q. Its oc- 
currence around 1 MHz is inconveniently close to 2.5 MHz 
in terms of filtering. In one scheme, this mode is suppressed 
by using two oppositely poled pieces of PZT to cause a 
zero net center-of-mass motion. This results in an order of 
magnitude reduction of the mode as measured by the real 
part of the input impedance. Shorter ring-down times, less 
cross-coupling, and more uniform acoustic radiation pat- 
terns result. 

The Lamb wave is suppressed by choosing a very thin 
foiL Since the asymmetric mode is largely responsible for 
Lamb- wave propagation, making the foil thin reduces the 
group velocity and prevents effective coupling from the 
array elements to the foih 

Since the Rayleigh wave is bound to within a few 
wavelengths of the backing surface, making deep cuts ex- 
tending well into the backing effectively suppresses it. 
Element Angular Response. The angular response of a com- 
pensated single element is shown in Fig. 7. This measure- 
ment was made by sending a signal to a point reflector at 
angle d from the normal to the transducer and then record- 
ing the return voltage. For such a round-trip or in/out type 
measurement, the normalized signal s[^] depends on fre- 
quency and element width as s(0) = 2Olog((sinx/x)cos^l 
where x = (ujL/2v] sinB, qj — radian frequency, and v = acous- 
tic velocity in water. The main points to note are that the 
element response drops off because of the finite element 
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width L and the obliquity factor cosfi. 
Lens, The acoustic beam generated by the ooe-djmeasionai 
phased array described here can oaly be electronically fo- 
cused and steered in the azimuth plane (see Fig. 3). To 
improve the resolutioD and sensitivity further, a cylindrica] 
acoustic lens^^-^^ is added to give a fixed locus in the ortho- 
gonal elevation plane. 

The main properties of interest in a !ens are the resolution 
at focus, radius of cur\^ature to achieve a given focus, the 
depth of field about that focus, and the intensitj^ gain 
caused by focusing. The lens must be sufficiently strongly 
focusing to achieve better resolution than an unfocused 
beam, yet have enough depth of field to allow imaging over 
a useful range. This involves a compromise in the elevation 
aperture W. In addition, reverberations caused by a poor 
lens impedance match to the patient are to be avoided. 
Finally, low attenuation in the lens and a convex outer 
surface are important considerations. 

The desired focal length F is related to the lens radius 
of curvatine R by R = F[l - l^nJ7i where yi = lA (empirically 
determined for the acoustic lens material used] , n = v^/v^, 
v^ - acoustic velocity In the iens, and Vj = acoustic velocity 
in the patient. 

The elevation be am width D^, at F is given by Dy = ^fFW] 
where \ = acoustic wavelength in the patient. The depth 
of focus D3^15Dy(l-0.01W/2F}. 

Array Fabrication 

Referring again to Fig. 3, the fabrication of an array probe 
involves the following processes, materials, or steps. 

■ Synthesis of a highly attenuative backing material with 
close impedance match to PZT 

■ Cutting, polishing, and plating the various ceramic ma- 
terials and matching layers 

■ Bonding the array together with sufficiently thin bonds 

■ Dicing the elements 

■ Providing electrical connections to the closely spaced 
elements 

■ Applying a top ground foil electrode 

■ Synthesizing and attaching a suitable lens 

1 Packaging the array into a probe complete with cable 

and connectors. 

The backing is required to match PZT acoustically and 
have sufficiently high acoustic loss to damp any back wave 
to a level below the system noise. Typically a 50'to-60-dB 
round-trip loss at any signal frequency above 1 MHz is 
required. If the loss is not large enough, multiple echos 
will return and cause bright lines or "range markers" in 



the image. The backing material chosen is a composite of 
tungsten for high impedance and polyvinyl chloride (PVC) 
plastic as a binder and dissipative medium. T>T}icaI param- 
eters achieved are an acoustic intpedance of 25 x lO^ rayls. 
a density of 14,2 x lO^ kg/m^. and an attenuation of 20 
dB/cm at 1 MHz. Typically the attenuation increases as 
frequency squared^ so at 2 MHz the loss is 80 dB/cm. 

The raw PZT-5H material used for each element is wa- 
fered to the appropriate dimensions, then ground and 
polished to a final thickness of 600 to 700 ^m. The thick- 
ness and the parallelism of the sides are carefully con* 
trolled. The parts are poled by a high electric field in a 
heated oil bath to make them piezoelectrically active. 

The backing and ceramic parts are very carefully cleaned 
under dust-free conditions, then bonded together with low- 
viscosity epoxy resin. Matching layers are also bonded at 
this stage» as is the electrical lead structure. To avoid the 
various cross-coupling problems outlined earlier, the 64 
indi%ddual elements are separated by sawing. 

The array is thoroughly cleaned again and a thin metal 
foil is bonded on to complete the ground connection and 
isolate the array elements from infiltration by grease, con- 
piing gel, or water. 

A variety of acoustic lenses can be used* A convex single- 
element lens is an effective choice^ but the selection of 
useful materials is limited to those in which the acoustic 
velocity is less than that in human tissue (1.54 x lo^ m/s). 
The lens material used for 77020A transducers is a 
polynre thane rubber, Sylgard 170®. It provides a good im- 
pedance match to human tissue and the acoustic velocity 
in Sylgard 170 is 1.02 x 10^ m/s. This allows the design of 
a convex lens with a radius R= - 36 mm for a 7-cm eleva- 
tion focal length. It is fabricated on a thin foil by casting 
and then Is bonded to the array. 

The final array is wired to a suitable coaxial cable con- 
taining 64 individual leads. It is then packaged in a water- 
proof plastic case. Fig. 1 on page 22 shows a photograph 
of a completed unit. 
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An Acoustic Transducer Array for Medical 
Imaging— Part II 



by David G. Miller 



HIGH SENSITFVrrY and a short pulse are important, 
but conflicting requirements for an acoustic trans- 
ducer used in diagnostic imaging systems. The de- 
sign of the transducer {Fig. 1) for the HP 77020A Ultrasound 
Imaging System is a compromise between these two re- 
quirements. Optimizing this compromise required careful 
consideration of the physical parameters and design con- 
straints discussed in the preceding article (page 17). 

Acoustic Pulse Propagation 

Lead zirconate titanate, a ferroelectric ceramic crystal, 




Ftg. 1 , One of the ultrasonic transducer probes used by the 
HP 77020 A Ultrasound imaging System, 



was chosen for the 64 piezoelectric elements in the trans- 
ducer array because of its high electromechamcal coupling 
and high dielectric constant. The piezoelectric behavior 
occurs as a result of either a zirconium or titanium ion 
locating in a stable off-center position inside the otherwise 
symmetrica] lead and oxygen unit crystsi cell [see Fig. 2). 
Because this ion has a large charge, this creates noncoinci- 
dental centers of negative and positive charge in the unit 
cell, forming an electric dipole. 

Polarizing the material orients a majority of these di poles 
in the same direction. This is done by applying an external 
electric field to crystals heated above their Curie tempera- 
ture and then cooling the crystals. This process, called 
poling, uses the property that above the Curie temperature, 
the unit cell is a body-centered cubic structure (Fig. 2a]. 
As the crystal cools, the ion in the center of the structure 
moves off center in a direction influenced by the applied 
electric field. After the crystal is cooled to room tempera- 
ture and the electric field is removed, the ions stay in their 
off-center locations to form the polar i/>ed crystal 

If a voltage is applied to opposite faces of a polarized 
crystal, a field is produced in the crystal that acts upon 
the dlipoles to cause expansion and contraction of the crys- 
tal. By applying an alternating voltage, acoustic waves are 
propagated from the stress discontinuities occurring at each 
face of the cn^stal as it expands and contracts. 

Each element of the 64-element transducer array is a 
0H5-mm-thick slab of lead zirconate titanate polarized in 
the plane of the slab and having electrodes on each end as 
shown in Fig. 2c on page 18, Compression and rarefaction 
longitudinal waves are propagated in both the forward and 
reverse directions. Waves traveling in the reverse direction 
are absorbed by the transducer's backing to prevent echoes 
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Rg, Z, The lead zfrconate titanate crystai unft cefl is a body- 
centered cubic structure (a) above the Curie temperature. 
Beiow the Curie temperature, the zircor\ium ion moves off 
center to form a tetragonal structure (by 

from objects befaind the transducer. This backing also helps 
damp the multiple reflections occurring inside each ele- 
ment. Waves radiating from such a transducer are rep- 
resented in Fig. 3. 

The acoustic impedance of the transducer elements is 
about 25 times greater than that of the human body. A 
sevenfold Increase in the magnilude of the acoustic echo 
from a target in the body can be achieved by reducing the 
impedance mismatch between the transducer elements and 
the body. This is done by placing an intermediate imped- 
ance-matching layer between the elements and the body. 
To prevent reverberation inside this layer and maximize 
the magnitude of the acoustic wave> the matching layer is 
made a quarter wavelength thick at the transducer's center 
frequency.* 

•Quarter-wavelerigiJi matching was fij^t used by early agtronomars who disj^ovefBd Ehai 
iatnl stars were brighter when viewed with ^eiescopss that had crfd tarnished lenses rarhef 
t^an riBW unlamished ones. "Hie lafriish layer lymsd oul to be about a quaJler wavelenglh 
thick ant] hsve an Int&rm^iate index ol rafractiorn betw^gn thut d the gl^as end the air 
In modem times, coaling optical coffiporkenta with qyarrsr-wave length-thick layers has 
iDecome commonplacB and interest in quarter- waveier>gth impedance malchingof etectrbal 
trensmi^ssions lines Tia5 resulted in an extensive developrr^nt oi the (tieory 
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Rg. 3. Longitudinaf compression and rarefaction waves 
traveling to the nght with velocity v for an acoustic transducer 
radiat!r}g pulses of wavelength k. Particles along the axis of 
propagation oscillate to the right and left. 



The constructive interference produced by the quarter- 
wav^elength matching layer is illustrated in Fig. 4. The 
acxiustic wave pul&e is represented as a ray r which strikes 
the layer at norraaJ incidence [shown at a slight angle to 
show the reflection r^ w^ithin the layer) and emerges as the 
sum of rays r and r^ . The magnitude and phase of r^ may 
be determined from the reflection coefficient 



R = 



^2-^ 

Z2+Z1 



(t) 



Thus a 180" phase change occurs for reflected waves 
when a pulse propagates from a high- impedance (Z^) region 
to a low-impedance (Z^J region. And it follows that since 
the reflected ray r^ travels an extra distance of one haU 
wavelength within the intermediate matching layer, this 
180° phase change aligns r^ to be in phase w^ith r so that 
the original pulse and reflected pulse amplitudes add as 
shown in Fig. 4. 

Equivalent Circuit Model 

To select the optimum acoustic impedance for the inter^ 
mediate matching layer and try to account for the effects 
of the electrode metallization and other features of the 
transducer array construction, an equivalent circuit model 
was developed. The KLM model ^ of an acoustic transducer 
is an electromechanical equivalent circuit (Fig. 5) in which 
the transducer is excited at the center. The acoustic load 
impedances can be transformed to the center node by the 
transmission line equation 



Ztn — Z( 



ZLCOsC-rrd/X) + jZ^ sinlirdlk] 



(21 



'^ Zp cos(iTd/X) + jZl sin('TTdyx) 
The circuit regards force (stress times area} as equivalent 
to voltage, and acoustic wave velocity as equivalent to cur- 
rent. 

Consider what happens when the front and back termi- 
nations are not perfect. The waves produced at the center 
suffer reflections at the faces, giving rise to a total stress 
and velocity that depend on the reflection coefficients at 
each surface. The stress, and hence the voltage V, will tend 
to be maximum when both the incident and reflected waves 
are strongs This occurs near resonance; therefore, the trans- 
ducer should be a half wavelength long for strongest exci- 
tation, if at resonance the crystal is a half wavelength long 
and each transmission line becomes a quarter wavelength 
long, Equation (2) simplifies to 



ZrvT — Zirt fAl 



(3) 



However, reflections from both surfaces tend to make the 
transducer ring, producing an output pulse like a decaying 
sinusoidal waveform whose decay rate depends on both 
the acoustic impedance mismatch and the electrical load- 
ing of the transducer. 

Denoting the transformed front and back load imped- 
ances by ZpQ and ZgQ, the circuit can be redrawn as shown 
in Fig. 6. Hence, the fraction of the power radiated to the 
body is 



Pf = 



V^/Zro _ Zr 

V^/Zbo + V^/Zpo Zb + Zp 



OCTOBER 1983 HEWLETT-PACKARD JO^JR^JAL23 



)Copr. 1949-1998 Hewlett-Packard Co. 




Unrefl^ct^ Fui$e r 



and the conflict between high sensiti^dt^'^ and a short pulse 
is apparent. High sensitivity requires a low backing imped- 
ance Z|j, but the shortest pulse requires matching the front 
and back acoustic loads to the crystal to reduce ringing. 

Matching the front of the crystal is further complicated 
because the crystal has an acoustic impedance Z^ of 
33x10^ rayls and a body's acoustic impedance is only 
1.5 X 10^ rayls. Thus^ the power percentage radiated to the 
body by a tiansducer with a matched backing and un- 
matched front is only 4%. 

The effects of the quarter- wavelength matching layer dis- 



Reflected Pulse r. 



Sum r+i-T 



Fig, 4, ff lustration of const fuctive 
interference in a quarter-wave- 

tength-thfok matching layer. The 
normal rays r and r^ are drawn at 
a slight angle to their true direction 
so they can both be seen. 



cussed earlier can be modeled by adding another transmis- 
sion line to the front port of the KLM model^ This layer 
acts as a mechanical transformer, causing a large load to 
be reflected to the transducer. With this matching layer, 
the radiated power percentage at resonance increases from 
4% to 27% and the same improvement occurs in receiving 
the acoustic echo. By using Equation (3), the acoustic im- 
pedance required for the matching layer can be calculated 
from 

2^IN = Zq = Z lay Ej^^ trod y 



Transducer Test System Design 



by Gsorge A. Fisher 



Early in ttie design of the phased-array transducer for the 
77020A. ttie need for accurate acoustic measurements was 
realized. As new designs were investigated, accurate and repeat- 
able measurements of the acoustic beam profile, main beam 
roll off, pulse length, ring down, and overafl array sensftivity were 
especiaHy important because these factors are critical to the 
overall performance of a phased- array imaging system. Measure- 
ments of individual element senstt^viiy and electrical impedance 
as a function of frequency provide additional information aboul 
acoustic beam performance, Individual element capacitance 
measurements provide a quick go/no-go check of array integrity 
for manufacturing and physical shock survival experiments. 

To perform this wide variety of measurements easily, a special 
test system is required. It must be flexible so that new experiments 
can be added easily as d ifferent areas of transducer performance 
are evaluated. It is also desired that the system be easy to use 
so that measurements can be made by lab technicians instead 
of engineers. Finally because of the huge range of signal levels 
resulting from echoes from various targets in the body, the total 
dynamic range of the system should be at least 80 dB 

To meet the above performance specifications, a test system 
(Fig. 1) was designed around an HP 9625A Computer With 
proper I/O expansion, this provides flexibility for new expenments 
as well as the intelligence required to provide ease of use and 
rrieasurement accuracy For mechanical accuracy, a special 
mechanical stage was designed to hold and position the array 
accurateiy in a tank of water along with a vanety of lest target 
reflectors and transducers. The sysiem also uses a regular pro- 
duction 77020A Ultrasound Imaging System scanner to control 
the firing and receiving of the array signals. A few custom 
hardware components were added to provide an interface to the 
scanner and tank positioning system, A flexible disc drive was 



added for program and data storage, and a plotter is used as 
the principal output device for analytical measurements. For im- 
pedance and capacitance measurements, a network analyzer 
and an LCR meter are interfaced to the 98 25 A 

The mechanical stage holds the transducer securely \n posi- 
tion, facing downward into a tank of water, which closely approx- 
imates the acoustic impedance of the body. The virtual central 
axis of the transducer is manually aligned with the axis of an 
adjustable arm that swings through an arc of ±50° under control 
of the 9825A The ann can be adjusted to position the target 
reflector at various depths of interest The stage holding the 
transducer is also controlled by the 9S25A and can be pivoted 
in the elevation direction. These movements allow positioning 
the target and transducer in azimuth and elevation directions for 
measuring beam directivity and focus. 

The 9825A Computer is the heart of the system. It provides 
convenient control of many variables during tests, as well as a 
simple operator interface. In addition, it enhances the accuracy 
of the system by performing calibration routines and then saving 
correction data. 

Another function performed by the 9825 A is that of extending 
the dynamic range of the system. The gain of the ultrasound 
scanner used to transmit and receive the acoustic beam is con- 
trolled by the 9825A by means of a digital -to- ana log converter 
(DAC), which feeds the TGC (time gain compensation) inputs of 
the scanner The gain of the scanner can be varied by 60 dB 
Before tests are performed, the 9825A feeds a lest signal through 
the scanner and measures the scanner's gain as a function of 
TGC voltage. This data is saved on a flexible disc During actual 
tests, the 9825A controls the TGC voltages to the scanner \n 
such a manner as to maintain the widely varying received signal 
in the middle of the available dynamic range of the scanner, and 
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Eiectrtcal Pott 

Rg. 5. Model of piszoefectric transducer after Knmhoiti, L ee- 
dom, and MatthaeL^ Here A is the cross-sectionaf area of ifie 
e/emeni. d is rfs iength, e^ is the clamf^ed dieiecjfic consmr}!. 
and Co = MJd is the damped capacitance C represents 
th^ acoustic eriergy storage, ifr/s the ratio of induced acoustic 
veiocity (equivalent current) to the efectricai current, Z^ is the 
characteristic acoustic impedance of the crystal, <i>(, is me 
haif-wavelength frequency of the crystat. Vf, and Vg ^^e stiff- 
ened acoustic veiocities. and Zf and I^ are the front and 
back acoustic toad impedances. 

which yields a layer impedance of 7.0 x 10® rayls. 

Lens Design 

The performance of the transducer is also determined 







FIfl, 6. Simplified edutvalent drcuu ustng transformed toad 
impedarKes. 

by the shape of the acoustic beam in both the direction of 
the scan and normal to the scan [elevation]. An acoustic 
lens was designed for the elevation focusing. 

The minimum angular beam width for a focal point 
placed anywhere in the range is given by the diffraction 
limit $ = KPN, where X is Ihe acoustic wavelength and W 
is the aperture.'* The geometric focal length of alens caji 
be derived by a geometric construction based upon Fer* 
mat's Principle* and the same approximation that produces 
the lens maker's equations for optical lenses.^ Thus, for a 
thin plane convex lens. 



Focal length = 



Radius of lens surface 



I * body' ^ lens I 



1 



' Formal 's Principle stales that a light ray wiEl take tfie patfi Detw.een tws points ihat fec|uife^ 
the leasl travel lime 



corrects the measurements by the stored gain values. This gives 
the system a total dynamic range considerably greater than 80 
dB 

For iTteasurements of both overaN and [ndivjdual element sen- 
sitivfty. the 9825A communicates with the scanner via the HP-IB 
{tEEE 488). putting the scanner into a special test mode that 
aJlows the com pule r ro take over control of transmit firing and 
receive phasing The computer fires just two elements in the 
center of the array, which results in nearly omnidirectional acous- 
tic radiation. This is reflected by a small target The computer 
also controls the scanner so that the signal from only one element 
at a time is processed and measured Thrs allows the serisftivity 
of each rndtvidual element to be measured. 



Peak RF 



The electric impedance and capacitance of each element are 
measured with the network analyzer and LCR meter Each instru- 
ment fs connected to one transducer element at a lime through 
a programmed relay matrix. This matrix contributes stray capaci- 
tance and inductance to the measurements, so before tests are 
run, the 982 5 A measures these stray values by performing tests 
on known impedances, calculating the errors, and storing these 
values For impedance measurements, the computer controls 
the measurement frequency with a DAC connected to the FM 
Input of the nelwork analyzer. At each frequency, and for each 
element in the array, the computer relates the measured imped- 
ance of the element to its actual impedance via a two-port network 
calculation using the stored data 




MP 9d25A 
Computer 



Fig. 1. Btock diagram of trans- 
ducer test system. 
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Fig, 7. Cafcutated round-trip beamwidths for 3S-MHz, fifxecf- 
foous, short-range and medium-range transducers. 

where acoustic velocities Vbo^y ^^^ ^uns ^^^ ^^^^ directly 
in place of the refractive indexes used in conventional 
optics. 

But geometric optics theory gives very poor beamwidth 
predictions because the apertures are relatively small — 10 
to 20 acoustic wavelengths. These small apertures produce 
beams dominated by diffraction effects and are not well 
described by converging and diverging rays. Another com- 
plication Is the pulsed operation rather than the continu- 
ous-wave operation for which the diffraction equation is 
valid. 

The optimum aperture and focal length combination pro- 
duces the mmimum beamwidth through the range of in- 
terest* Two different combinations are illustrated in Fig. 7. 
The short-range transducer has a short focal length lens 
and a small aperture to reduce the beamwidth near the 
transducer. The medium-range transducer uses a larger 
aperture to reduce the beamwidth at midrange. 

The acoustic lens design also requires the selection of a 
material with the proper acoustic velocity and impedance 
to match the human body. Velocities lower than that in 



the body allow a convex face on a flat transducer array. 
Velocities higher than in the body require a concave lens 
surface or the use of two lens elements. Other lens require- 
ments are low attenuation for all beam steering angles, 
resistance to abrasion, chemical resistance to acoustic cou- 
pling gels, cleaning disinfectants, and sterilization, and 
providing a permanent seal to the transducer package. A 
convex lens is used for the 77020A's transducers* 
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Radiated Power Characteristics of 
Diagnostic Ultrasound Transducers 



by Thomas L. Szabo and Gary A. Seavey 



WHAT LEVELS OF ULTRASONIC POWER are 
safe? This is an often asked question to which 
there is still no absolute answer, but some 
guidelines have been developed. People working with diag- 
nostic imaging and other uses of ultrasound in medicine 
have established safe levels of ultrasonic intensity below 
\uhich no bioeffects have been demonstrated. Ulrich^ de- 
veloped a threshold based on a compilation of biological 
effects data. He arrived at a crilenon for minimal hazard: 
a spatial-average, temporal-average intensity below 100 



mW/cm^. This criterion and others have been used to sum- 
marize the lowest levels at which bioeffects have been ob- 
served. For example, see Fig. 1 adapted from Nyborg-^ All 
of these observations can be sunnnarized by the following 
statements: 

■ No substantial bioeffects have been demonstrated for 
spatial-peak, temp oral -average intensities <100 mW/ 
cm^. 

■ No substantial bioeffects have been demonstrated for 
exposures for which the product (Ixt) is less than 50 
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Fig* 1i Suggested safe uitrasound mtenstty ieveis based on 
the kmes! ieveis at which bioiogfcal effects have been re- 
ported fn tissues (adapted from reference 2). 

J/cm^p where for pulsed operation, t is total on-and-o£f 

time, and I is the spatial -peak* temp oral -average mten* 

sity. 

These results are incorporated in the "Safety Standard 
for DiagnosUc Ultrasound Equipment/*^ Definitions and 
measurement recommendations outlined in the standard 
are used by the HP Medical Group for characterization of 
ultrasound transducers and systems. Before explaining 
these measurement techniques, it is necessary to answer 
the question: "How does a transducer provide a certain 
level of acoustic intensity in the body?" 

Transducers 

A transducer serves two main functions: converting elec- 
trical drive signals into acoustic waves, and directing and 
focusing these waves in the body. The transducer also re- 
ceives acoustic echoes from the body that are used to form 
the ultrasound image, but the emphasis here will be on the 
transducer as a transmitter. 

Although the HP 21205B 3.5-MHz Transducer (see Fig. 
1 on page 22] appears externally as one transducer, it is 
actually a phased array of 64 transducer elements. Each 
element has a certain electro acoustic conversion trans - 
ducer efficiency TE* when connected to a continuous-wave 
(CW) voltage generator V^^ with an impedance Rg. In terms 
of acoustic peak power Pg at the face of the element, TE' = 
P„/Pg, where P^ ~ Vg^/4Rg is the generator power. 

A transducer can be thought of as a bandpass filter with 
a certain impulse response. To find the temporal-average 
acoustic outputt the input waveform also must be consid- 
ered. Most diagnostic imaging transducers are driven by 
short pulses repeated at long time intervals. A common 
drive pulse is a tone burst of n half cycles* each t^ in length. 
The shape of the input pulse determines the temporal- 
average power per pulse. For a tone burst, this power is P„/2. 
The average power is half as large as the peak power, so 
an average power factor apf can he defined in this case as 
apf = 0.5. Another important factor in reducing total aver- 
age power is the duty cycie, which is the drive pulse dura- 
tion Ip divided by the pulse repetition interval Tj^, In the 
HP 77020A Ultrasound Imaging System, the duty cycle is 
typically 0.0005. Fig. 2 compares a computer simulation 




Lffi 13 I4e IM e.88 LBB L2B Hi \M LBB IM 

Time {$ib} 

Fig, 2. Measured and computer simulated temporal pres- 
sure waveforms at the spatiai peak for the HP 21205B Trans- 
ducer. 

of one acoustic temporal pressure output waveform of a 
single 3,5-MHz element in water with an actual waveform. 

The temporal -average acoustic power of all of the ele- 
ments at their faces can be calculated using the equation: 



P_ = N(TE)(P„apf](tprrR) 



(1] 



where TE is now defined as the temp oral -average acoustic 
output divided by the average electrical input power for 
the drive pulse used. 

Beamformfng 

Tu improve image definition, the trans ducer array directs 
and focuses acoustic energy by two methods. With refer- 
ence to Fig. 3, the beamformingcan be separated into ortho- 
gonal elevation and azimuth planes. In the elevation plane 
perpendicular to the transducer element widths, a fixed 
lens covers the array and focuses the beam. At the focal 
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Fig. 3. Definttion of acoustto beam profiies in azfmuth and 
etevaUon pfanes. 
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plane Z = f , the beam shape for an elevation aperture L^ is 
the Fourier transform of the constant-amplitude aperture 
hinction sin ('7TLeX/FX)/('TrLeX7'FXJ. A simple estimate of this 
large compression of the beam is the haif-amplitude (6-dB) 
half heamwidth. 



Ti'anaaucef 



x^jB = l,9FK/TrU 



(2) 



where X is the acoustic wavelength in water. 

In the azimuth plane, focusing and beam steering are 
accomplished by addressing each phased-array element 
with appropriately phased, time-delayed signals. As shown 
in Fig. 3, the beam is focused in azimuth along the Z axis 
at the same place as the elevation focus. The Z axis consti- 
tutes one line in the ultrasound image (0° line); other lines 
are formed by directing the beam to other angles in the 
azimuth plane* At the focal plane, a 6-dB or quarter-power 
ellipse can be defined with half beamwidths (Equation (2j) 
for semiaxes. If the spatial -average, temporal -average inten- 
sity Is ATA is defined at the focal plane as the average power 
in the integrated beam cross section divided by the area of 
this ellipse, the following formula applies to the ideal 
example discussed. 



1. 



°Fgvg^titil^a 



m 



where L^ is the total azimuth aperture. When the intensity 
at the peak is measured by a microprobe transducer and 
time averaged, the value is referred to as the spatial-peak, 
temporal-average intensity Ispta- 

Another commonly quoted intensity is the spatial-peak, 
pulse-average intensity Isppa. This single-pulse intensity 
is the temporal average without the duty cycle t^/Tj^. There- 
fore, IgppA = IsptaC^h^^p) ^s evident from Equation (i). 

To compare actual measurements with predicted values, 
a measurement scheme, outlined in reference 3, is used. 
This method uses a hydrophone probe, its calibration, and 
integration algorithms for spatial and temporal averaging. 

Transducer Hydrophone 

To characterize the radiated acoustic power both spa- 
tially and temporally in absolute power terms, a special 
calibrated miniature hydrophone is used. This hydrophone 
has a piezoelectric element approximately 1 to 2 mm in 
diameter fabricated of polyvinyiidene fluoride (P VF2) film. 

Since acoustic intensity is proportional to the pressure 
squared, and the hydrophone output voltage V is propor- 
tional to pressure, a hydrophone sensitivity can be defined 
as K = I/V^, where 1 is the instantaneous intensity and V 
is the instantaneous voltage response of the hydrophone. 

An X-cut quartz transducer calibrated by the U.S. Na- 
tional Bureau of Standards is used as the reference power 
source to determine the hydrophone's sensitivity. The 
radiated acoustic power P^ versus transducer drive voltage 
Vg has been measured using a force balance to arrive at a 
radiation conductance Gj., where Pt = GrVg. The hydro- 
phone is positioned in the far field on the beam axis of the 
quartz transducer* The hydrophone's response V(0,t] to the 
radiated temporal pressure is recorded. A number of radial 
scans are taken to characterize the spatial distribution of 
the radiated pressure. Since the total radiated acoustic 
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Water Bath 
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Fig. 4. Acoustic intensity measurement apparatus. 

power from the quartz transducer is known from P^ =GrVg, 
knowing the temporal and spatial characteristics of the 
radiated pressure, one can do a temporal and spatial aver- 
age of the pressure squared to relate the hydrophone's sen- 
sitivity to the measmement taken. This intensity response 
factor K for a given frequency is 



K = 



hb 



m 



where Pj is the total radiated acoustic power in watts from 
the source transducer into the water, Vp(0) is the temporal- 
peak hydrophone response on the beam axis, and h is the 
temporal average of the hydrophone response on the beam 
axis. That is, 



Tr Jo 



(o,t]dt 



(5) 



where V[0,t) is the hydrophone response waveform. Here 
b is the spatial normalization computed from the spatial 
scan measurements* For rectangular transducers. 






Vp(x,0]VptO,y)dxdy 



(6) 



where Vp(x,0] and Vp(O.y) are temporal-peak hydrophone 
responses along orthogonal beam profiles through the beam 
axis, and x^^q and yg^^Q are the 6-dB half beamwidths. 
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Fig* 5. Acoustic pressure profits or beamplol at the focai 
piane of 3 transducer compared to a probe mBasurement 
simulation. 
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Measurements and Deflnitfons 

To clmracierize the radiated acoustic power from an ul- 
trasotLEid imaging system, one must be able to find the most 
intense part of the transducer beam m three dimensions 
and measure the acoustic pressure waveform there, A 
sketch of the water tank arrangement used to make these 
measurements is shown in Fig. 4. The support structure 
(not shown) allows accurate aligoment of the transducer 
and hydrophone and automated positioning of the cali- 
brated hydrophone. 

When the hydrophone is positioned on the beam axis at 
the spatial peak, the calibrated hydrophone temporal out- 
put waveform V(0,t) is recorded. A t\T]ical waveform is 
shown in Fig. 2 and it agrees well in shape with the pre- 
dicted pressure waveform for a single element. The mea- 
sured waveform is squared and integrated over a long 
period Ti^to determine Ispta using the following equation. 



I. 






t)dt 



(7) 



To find the spatial average of the intensity, more detailed 
information about the beam is gathered. To do this, only 
the temporal value Vp is used from each waveform for these 
spatial measurements. Next, the hydrophone is scaiuied 
along two orthogonal paths within the spatial- peak plane 
to map the pressure profiles, Vp(x,0] and Vp(0,y). 

A typical spatial pressure profile such as Vp(0,y) appears 
similar to the ideal pressure beam profile in Fig. 5, At the 
6*dB points in each pressure profile, a half beamwidth is 
determined. These half beam widths, x^j^g and y^^B' ^^ke 
the ideal of Equation (2), determine the measured focal 
ellipse depicted in Fig. 3. This ellipse has an area A = 
TrXg^jByedB- The final value of Isata is given by 



U 



SFTA 
A 



V^(0,0] 



(8J 



This equation is based on measurable values and its result 
can be compared to the ideal estimate given by Equation (3). 



Conclusions 

At this time, safe levels for diagnostic ultrasound systems 
are those agreed on by consensus. Carefully worded defini- 
tions insure that those who make intensity measurements 
are speaking the same language. Results for diMerent tyi^es 
of diagnostic systems are listed in Table 1. HP systems fall 
in the range of an automatic sector scanner in this table. 

HP transducer arrays and ultrasound systems are 
evaluated according to the methods described above. 
Transducer modeling and beamplot computer programs, 
along with equations such as Equations (1) and (3), have 
been used to predict acoustic intensity^ levels with good 
accuracy. Hence, the mechanisms that produce acoustic 
intensities are well understood and can be changed by 
design. 

Work is also continuing on refining the measurement 
process itself. In Fig. 5 the effective hydrophone probe 
diameter is showm in relation to a typical pressure profile. 
Because the beam is comparable in size to the probe diam- 
eter, the probe smooths or integrates the actual pressure 
profile. The estimated hydrophone response is the second 
curve in Fig. 5. Because the hydrophone probe was cali- 
brated originally against a quartz source much larger than 
its diameter, its calibration factor K is strictly applicable 
only in similar environments. The effect shown in Fig. 5 
is that a slight error of about 1 dB is introduced at the 
spatial peak. Other sources of error are under investigation 
to improve the measurement accuracy further- 
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Table I 
Intensity Ranges Produced by Current Ultrasound Syetems 



Type of Equ ipment 



Spatiel-AveragB, 
Tern poraj -Average 

[ntBn&ityontbe 
Radiating Surface 



Spatial'Peak, Temporal- 
Average Intensity 



Spatial'Peak, PuIb&- 
Average Intensity 



Static PylseEcho Scanners anil M-Mode Equipment 0.4-Za mW/cro^ 

Automatic Sector Scanners ( Phased Arrays and Wobblers) Ji.7-80 mW/cin' 

Sequenced Linear Arrays 0^06-10 mW/con^ 
Pulaed Do p pier. Primarily for Cardiac Work 3-32 mW/cro^ 

Doppler Instmments, Primarily for Obstetric Applications 3-25 mW/cm^ 

CW Doppler, Primaf iiy for Peripheral Vascular Applications 3B-84D mW/cra^ 



10-200 mW/cm* 

45-200 mW/cm^'* 

0.1-UmW/cm^ 

50-290 mW/cm^ 



0.5-280 W/cm* 

25-100 W/cm^ 

25-100 W/cm^ 

3-14 W/cm^ 



Spatial Peaklntensityi 9-75 mW/cm^ 
Spatial Peak intensity: 1 1 0-2500 raW/cm'* 



Estimale based on the spstlal average — temporal average* 

This value was measured with the scanning mechaniam arrested for M-mode end at the maximum system putse repetition rate. 
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A Scan Conversion Algorithm for 
Displaying Ultrasound Images 

by Steven C. Leavitt, Barry F. Hunt, and Hugh G. Larsen 



IN THE HP 77020A Ultrasound Imaging System, the sig- 
nals received by the 64 acoustic transducer elements 
and processed by the electronic scanner conespond to 
the return echoes hom an acoustic beam focused in a direc- 
tion that changes with time in a predetermined manner 
The digital scan converter in this system processes this 
information to present a visual image of the acoustically 
scanned region. This converter uses a novel two-dimen- 
sional interpolation algorithm in a process referred to as 
R-Theta reconstruction. Using this technique^ the acoustic 
image is reproduced from the echo values in an essentially 
error-free manner. 

Review of Acoustic Data Processing 

After the 64 individual signals received by the transducer 
array are appropriately phased, delayed, and summed, the 
result is a single waveform consisting of bursts of sinusoidal 
activity- This activity is a fairly broadband set of pulses at 
frequencies near the transducer's resonant frequency. Spa- 
tiallyt this waveform corresponds to the echoes along a 
narrow^ line directed radially outward from the center of 
the transducer's face, A total of 121 such acoustic rays are 
used to form a 90° sector image (Fig. 1). The waveform for 
each ray is routed to a full- wave rectifier ^ a peak detector, 
and a low -pass filter to extract the modulation envelope. 
The envelope is sampled, first at 15 MHz. and then at a 
slower rate dependent on the selected scan depth. This 
allows the use of a fixed number (396) of sample points 
per ray to represent different scan depths of 4, 8, 12, 16^ 
20. and 24 cm. The array of 121 x 396 five-bit samples is 
sent from the 7 7020 A 's scanner to the digital scan converter 
for image formatting. 

Hence, the data received by the digital scan converter is 
the result of a two-dimensional sampling process. Not only 
is the data sampled radially, it is sampled laterally in the 
angular or azimuth dimension. All of the concepts relevant 
to the Shannon sampling theorem associated with a time- 
based waveform (the radial sampling) apply equally well 
to spatial sampling. The implication is that, given samples 
from a tw^o-dimensional sampling process, and if the sam- 
pling frequencies are sufficiently high, then in principle, 
the image may be reconstructed w^ithout error from these 
samples. 

If the azimuth sampling is viewed as a spoke density, 
then in this case, the spokes or acoustic rays are 0.75* apart. 
Small spatial events that consist primarily of high spatial 
frequencies relative to the spoke density will be missed. 
Furthermore, since the acoustic field is viewed through a 
finite window (the transducer) there is a resolution limit 
related to the transducer aperture and the acoustic 



wavelength. For this reason, it makes no sense to increase 
the spoke density if the aperture-related resolution does 
not justify it. 

Scan Conversion Methods 

Digital scan conversion is a relatively recent innovation 
in ultrasound imaging. Earlier systems used direct analog 
converters to modulate the intensity of an X-Y display. 
There are problems with this technique that outweigh its 
simplicity. First, a completely analog conversion results in 
a picture that looks like spokes of a wheel with an image 
superimposed on it. By today's standards, this is unaccept- 
able image quality. Also, because storage of the data in a 
digital memory allows high-quality freeze- frame display 
and the ability to manipulate the image using various gray- 
level mapping curves and digital enhancements, analog 
conversion is no longer a viable option for medical ul- 
trasound imaging. 

Many digital scan conversion methods have been de- 
scribed^ that map polar coordinate data to a rectangular 
format, but most of these suffer digital artifacts caused by 
the scan conversion process. To review some of these 
methods, the term scan conversion is used here to refer to 
the spatial transformation process by which data that natur- 
ally occurs in a polar coordinate system is mapped into a rec- 
tangular memory or grid. The problem is that the precise loca* 



Display Memory 



396 Data Samptas 
Atong EsctT Line 




See Fig. 3 



Sector Display 



Fig* 1. The display of a 9Cr acoustic SBCtor scan /s ac- 
comptished by mapping the 121 scan fines, each contafmng 
396 data samples, onto the HP 77020A Uftfasound Imaging 
System's B40 x 4B0-pixel display memory as shown. A mag- 
nified view of the enarded area is shown sn Fig. 3 to iitustrate 
the mapping problem. 
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tioD of each sample point in space does not, in general, 
correspond to an allowable display point. Therefore, the 

data must be adjusted or scan converted before it caa be 
written to the display memory. 

The typical digital scan converter maps the scanner data 
to the nearest display memory location, that is, the polar 
coordinates are converted to rectangular coordinates, 
which in turn are quantized to the nearest pixel {di5play 
point). Each ray or scan line maps to the memory' grid in 
a somewhat irregular fashion. For example, portions of the 
data from neighboring lines can overwrite each other. The 
resulting image appears as a series of straight line segments. 
This basically incorrect process* where the data is pushed 
and pulled into alignment with the memory grid, results 
in an annoying artifact called a Moire pattern, which is a 
well-defined pattern of holes in the image corresponding 
to the unaddressed pixels (Fig. Za). This pattern is often 
reduced by creating interpolated scan lines between the 
actual scan lines and converting these in an identical man- 
ner. However, the overwriting problem then becomes more 
severe. Another method fills the holes by horizontally 
smearing immediately adjacent data into each hole to com- 
plete the scan conversion [Fig. 2b), but this introduces 
other artifacts. Hence, most of the current scan conversion 
methods create an image with a "steppy and blocky" 
appearance, 

R-Theta Scan Gonversfon 

The 7 7 020 As R-Theta scan conversion method uses a 
real-time reconstruction technique. The intensity values 
for ail pixels that lie between two adjacent scan lines are 
interpolated using Equations (1], [2), and (3). 



2ro ~ 2q( ji J ( 1 Rkrr ) + ^Qf j^ 4 1 |Rerr ( 2 ] 

2 = 2rq( 1 - B^uk) + ^K\ ^ERR ( 3 ] 

These equations result from the evaluation of the 
generalized interpolation formula 



Z{r,0) = X 2 S[r- jAr, O-kAB]Z{]ATMA0) 



(4! 



where the interpolation function S(r.O) takes on the spe- 
cific form of a two-dimensional triangular function. 
ZHAtXAO] corresponds to an actual sample value and 
Z{t,B] is an interpolated value at the pixel site. Equations 
(Ih (2). and (3) are best explained by referring to Fig. 3. 
Zqj^] and Zj^ni are equivalent samples (in a radial sense) 
on two adjacent scan lines with Zij^j representing the sam- 
ple on the line with a shallower angle from the normal to 
the transducer. 

Equation (1] radially interpolates a new value Z^j be* 
tween the inner and outer sample points on the inner radial 
line based on the radial position of the selected display 
point (Z in Fig. 3). Equation (2) performs the same interpo- 
lation on the outer radial line (Zrq). Equation (3) linearly 




(e) 



Fig. 2. Compamon of displays resutUng from (^} a commonly 
used scan conversion algorithm, (b) horizontal fill-m enhanae- 
mer^fs to the process used for (a), and (c) the R- Thela conver- 
sion algonthm used in the 7702QA. The left half shows half 
of a concentric ring image and the right haif shows half of a 
radial spoke image- Note the Moir6 pattern present In (a) and 
the stepped tines in (b). 
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Outer Scan Ltne 



Display 
Memory 




Inner Scan Une 



Radial Arcs 

at Sample 

Sites 



Calculated 
"^ Radius of Z 



Calculated 
Angle of Z 



Fig. 3. Enlarged view of a portion of one slice of the sector 
Trie dots represent data samples received from trie scanner 
erjd ttie crosses represent the display pixels the data samples 
are to be rnapped to. See text for additional explanation. 

interpolates Zri and Zj^o i^ *h€ azimuth (O] dimension in 
a similar manner to assign the proper value to the display 
point at Z. The 77020 A's digital scan converter uses this 
interpolation algorithm to assign values to all of the display 
points bounded by two successive scan lines. The process 
is repeated for all scan lines in the displayed sector. 

Hardware Design 

The block diagram in Fig. 4 details the circuitry involved 
in implementing the R-Theta scan conversion. 
Input Buffers and Multiplexers, Since each radial scan line 
is received from the 77020A's scanner in an asynchronous 
fashion, the sync/iatch function is nothing more than a 
digital resyncKronizatiort. Each data or header word is ac- 
companied by a strobe which indicates when that word is 
valid. The word rate, and hence the strobe frequency, varies 
with imaging depth, but never exceeds 3,8 MHz. The strobe 
is sampled at 12 MHz, and when found true, the remaining 
eleven data lines are latched and then interpreted. 

One of the lines, the format bit, indicates whether the 
word is a header (line number information] or a data word. 
If it is a header word, a new^ scan line is imminent. Thus, 
the digital scan converter must be finished with the line 
it was processing. Also, the remaining bits of the header 
word indicating the line number [angle of the scan] must 
be placed into the 3-header FIFO (first in, first out) buffer 
for future use. In addition to indicating the start of a new 
line, the header word indicates the end of the previous 
line, and therefore causes the rotating 3-line data buffer to 
rotate, putting the currently filling buffer into read mode 
and preparing Lhe buffer with the ''oldest*' data to receive 



data. Receiving a header word causes a line number to pop 
out of the 3-header FIFO buffer. This line is processed next 
by the digital scan converter. The R-Theta error calculators 
are reset to zero and the X and Y address counters are 
preset to a processor-selectable '^origin." Finally, the X-Y 
raster state machine is reset to the apex of the scan slice 
and the new line number from the 3-header FIFO buffer ts 
passed to the R-Theta error calculator and the X-Y raster 
state machine, which need angular information to perform 
their functions. 

The rotating 3-Hne data and 3-header FIFO buffers work 
in concert so that, as a new line number is emerging from 
the FIFO, the data for that scan Une is also emerging from 
the rotating buffers. These provide the three levels of delay 
necessary to interpolate between two scan lines, because 
the header word cannot be processed immediately when 
it is received, but rather must be delayed for two line 
periods while the data is gathered. The 3;1 multiplexers 
are used to .select the inner and outer lines for input to the 
8 -stage shift registers. 

The dual 8 -stage shifters are essentially a cache memory, 
retaining the last eight data samples from each of the two 
hues for which the slice is being interpolated. These 
bidirectional shift registers allow sequential recall of previ- 
ous data samples* The first two stages of each shifter pro- 
vide the inner sample (closest radially to the sector origin] 
and the outer sample (more distant) for each of the lines» 
Thus, the outputs of the shifters are the four sample values 
that surround the display pixel for which an interpolated 
value is sought- As the slice rastering progresses from pixel 
to pixel, two signals cause the shifters either to receive the 
next sample from the buffers (INCR) or to back up to retrieve 
a previously used sample (DECR). Of course, the buffers 
must also back up when DECR is true or advance when 
INCR occurs. It may not be obvious why it is necessary 
sometimes to retrieve old data samples. Referring to Fig. 
3, which represents a slice near the 45"" line on the left side 
of the sector, it can be seen that as the raster path traverses 
from pixel A to pixel B, and then to C, an integral radial 
boundary (represented by the arc from Zo^n) to Zjjnj) is 
crossed. Hence, the fom samples surrounding pixel C axe 
not the same as those surrounding pixels A and B. The 
same situation holds true for the right side of the sector 
and for all scan lines with steep angles. 

To handle this situation, each time an integral radial 
boundary is crossed, an INCR or DECR is generated, depend- 
ing on the direction of the crossing. If the raster crosses in 
a direction that decreases the radius of the slice, for exam- 
ple from pixel B to pixel C in Fig. 3, DECR is true and the 
shifters are left shifted, retrieving samples Zyj^-ri' ^un]* 
Zo(n-i)' ^^^ ^0[n] to provide the new correct corner values. 
Samples Zjjn+n and Z^^^^ ^ are purged from the shifters 
by being shifted out the left end into a bit bucket. Howeverj 
the rotating buffers have also backed up so that ^M^n+i] ^^^ 
^0[n^ii ^re now^ present on their outputs. 

When the raster scan passes from pixel C to pixel D, the 
radial boundar>^ is crossed in the other direction, yielding 
an INCR. The shifters shift right, dropping Zm^-^] and 
2^ofn i) o^* ^6 right end, and the rotating buffers now 
make Z|,n + 2) ^^^ ^q^^^s] available at their outputs. Out of 
the first two shifter stages now come Zj^j, ^ ^^ Zii^j, Zq^h * iji 
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Fig. 4. Block diagram of the R-Theta scan diversion system. 



and ZQ[n|t which represent the four samples surrounding 
pixel D (as well as pixels A and B earlier). 
X"Y Raster State Machine. The interpolation starts with a 
large state machine programmed to scan from the apex of 
the slice down to the slice's outer radial limit in a predeter- 
mined way. Each slice has a unique traversal path selected 
to minimize addressing pixels not lying in the currently 
addressed slice and to simplify encoding. EPROMs store 
the encoded traversal paths and the paths used for the left 
half of the sector are reflected by symmetry into its right 
half. 

By using an extremely simple encoding algorithm, each 
EPROM can store eight slice traversal paths. With the slice 
traversal always beginning at the apex, a logic one indicates 
a drop In the Y direction and a logic zero indicates move- 
ment in either X direction. Note that only Y increments 
are used; X addressing is made bidirectional by using an 
additional flag bit XDIR, which toggles every time a move- 
ment in Y is made. The XDIR flag is set to one at the apex 
to indicate that the next X movement is outward to the 
outer scan line. This holds true for slices on both sides of 
the sector because of the definitions used for inner and 
outer scan lines. 

A STARTLN signal generated by a header word from the 
scanner initialises the state address counter. The output 
of this counter and the line number form the address for 
the EPROMs. The traversal bytes are read by the state ad- 
dress counter indexing through the 256-byte address space 



for that particular line number. The three least -significant 
bits of the state address counter serialize the byte data to 
form the YINC signal. Asserting YINC causes the memory 
address control to increment the Y address counter with a 
YCNT signal. Deasserting YINC generates^ XCNT signal, 
which is coupled with XUIR to form XUP/DN, which incre- 
ments or decrements the X address counter. 

The X and Y address counters are loaded after every 
STARTLN signal with the sector origin as specified by the 
microprocessor* This positions the sector data within the 
digital scan converter's B40X-480 memory space since all 
pixel movements dictated by the state memory are relative 
to the apex. Although this accounts for how the address 
counters sequence in a raster-like manner down through 
each slice, the interpolation pipeline still must be told in 
some manner which samples are to be used for interpolat- 
ing and what the interpolation weights are. 

The error weightings Rkrr ^^i"^ ^ehh ^^^ necessary to use 
interpolation equations [1], [2], and [3], Since the radius 
and angle of each sample is known from its header word 
and sample number, the straightforward approach is to 
calculate the polar coordinates of each pixel site within 
the sector. Then, with a little comparison circuitry between 
the piixei polar coordinates and the sample polar coordi- 
nates, the values for R^rr ^^^ *^rr could be calculated. 
Using the apex as the reference point for the rectangular 
coordinates, the equations that must be solved for every 
pixel addressed by the state memory are 
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(5) 
(6) 



R-Theta Error Calculalor. Because a new pixel is addressed 
every 160 ns. simplifications were implemented in 
hardware to calculate pixel polar coordinates Rp and Op, 
These caJculations are performed by the R-Theta error cal- 
culator. The inputs are the line number, YINC, and XDIR. 
The outputs are Rerr, O^rr, VAL, INCR, and DECK. Re^r 
and ©ERR 3Je fed to the interpolation pipeline for calculat- 
ing the interpolated data value. They represent the radial 
and angular distances between the currently addressed 
pixel and the innermost sample on the outer scan line. 
They are shown in Fig. 3 where Zoj^j is the reference. VAL 
is used to qualify the pixel as being within the slice bounds, 
that is, its radius and angle information satisfy the currently 
addressed slice. 

Interpolation PipeUne. Once Rerr and O^nf^ are determined 
and the four values, Z|(n)i Z,„, + ij, ^o{n]^ *ind Zofn + i]t are 
also available, the interpolation takes place. The interpola- 
tion pipeline, although requiring significant hardware, is 
a very straightforward process. The four corner values are 



radially interpolated to yield two intermediate values> Zrj 
and Zrq, according to Equations [1] and (2). Then the final 
value Z corresponding to the gray level assigned to the 
associated display pixel is determined by Equation (3). 

The final step is to write Z to the correct address in the 
display memory. This operation is regulated by VAL, which 
as explained earlier, is true only when the pixel is w^ithin 
the slice. If the pixel is not in the slice, but was in the slice 
traversal path, VAL is false, disallowing the final write to 
the display memory. 

The scan conversion system's pipelined architecture in- 
terpolates a new pixel every clock cycle (160 ns}. A typical 
slice is composed of 1026 pixels, and a 90° sector contains 
about 123,000 pixels. Each of these is interpolated once 
ever>^ video frame (16,7 ms]. Fig. 2 illustrates the improve- 
ment in image quality gained by using the R-Theta al- 
gorithm. 

Reference 

1 , J , Ophir and N.F. Maklad, "Digital Scan Converters in Diagnostic 
Ultrasound Imaging." Proceedings of the IEEE, Vol, 67. no, 4. April 
1979, 



Ultrasound Image Quality 



by Richard A. Snyder and Richard J. Conrad 



IMAGE QUALITY is a measure of the diagnostic utility 
of an ultrasound image to a physician. Generally the 
physician uses ultrasound to examine the internal or- 
gans of a patient noninvasively. Hence, the ultrasound 
image must be an accurate representation of structures in 
the body. 

The human body Is an ultrasonic medium that is far from 
ideal in its ability to propagate an undistorted ultrasound 
wave, because the body attenuates, scatters, and refracts 
the ultrasound wavefront. The factors influencing the qual- 
ity of the ultrasound image include resolution, depth of 
penetration, and tissue representation* 

Resolution 

The resolution of the ultrasound image is a measure of 
its ability to separate closely spaced individual structures- 
Referring to Fig. 1 . the resolution for a phased-array sector 
scanner is divided into three distinct areas: the resolution 
in the direction of the scan angle, or azimuth resolution, 
the resolution in the direction perpendicular to the plane 
of the scan, or elevation resolution, and the resolution in 
the direction of pulse propagation, or range resolution. 

The sector image is formed by transmitting an ultrasound 
pulse of a particular center frequency in one direction into 



the body and receiving and processing its echoes. The di- 
rection of propagation defines the position of one scan line 
of the image while the arrival times of the echoes determine 
where along that line the echo sources lie. This transmit/re- 
ceive process is repeated in different directions until the 
entire plane of the sector has been scanned. 

The intensity distribution in space of the transmitted 
pulse as it propagates is known as the transmit beam pat- 
tern. The sensitivity distribution in space of the receiver 
is known as the receive beam pattern. The total round-trip 
azimuth resolution and elevation resolution are each deter- 
mined by the product of these two beam patterns. The 
distinction between elevation and azimuth resolution is 
necessary because the aperture of the acoustic transducer 
IS rectangular, with one dimension falling along the eleva- 
tion direction and the other along the azimuth direction. 
These dimensions are not necessarily equal in length, so 
the round-trip beam patterns in the azimuth and elevation 
directions may be different. Each beam pattern is a function 
of the aperture size, the wavelength of the ultrasound pulse, 
and the focal distance of the acoustic lens. Because the 
aperture size is typically only 32 times the wavelength, or 
less, diffraction effects play a significant role in the analysis 
of the beam pattern. The first-order diffraction effect causes 
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Fig, 1 . An ultrasound image for a 9(f sector is ohtatned by 

scanning a beam of uftrasound fiulses m an azimuth direction. 
The beam's azimuth, elevation, and radml resolutions are 
functions of many factors as discussed in the text. 



the beam to diverge in the transducer's far field. The far 
field region begins at a distance Z = W^/X for a rectangular 
aperture, where W is the aperture size and X is the 
wavelength. In the near field, the beam width may be re- 
duced by focusing with a lens or focusing by controlling 
the phase of each element in the phased array. The beam- 
width at the focal point may be reduced and the near field 
extended by increasing the aperture. The price of this im- 
provement in resolution at the focal point is that the beam- 
width before or beyond the focal point increases rapidJy 
with aperture causing a shorter depth of fieldn [n fact, al- 
though beamwidth at the focal point decreases linearly 
with aperture size, depth of field decreases with the square 
of the aperture size. This same loss in depth of field is 
experienced when the focal point is brought closer to the 
aperture. This situation forces a compromise between res- 
olution at the focal point and depth of field whenever a 
fixed focus is used. Fig, 2 illustrates the approximate trans- 
mit beam patterns from transducers with apertures of equal 
size, but with four different focal points. Although the 
minimum width of each beam pattern does decrease as the 
focal point is moved toward the transducer, the beamwidth 
diverges more rapidly away from the focal point- 

For the phased-array scanner, the transmit focal point in 
the azimuth direction is controlled by the delay applied 
to each element of the array when the pulse is transmitted. 
This focal point is fixed at a depth that provides the best 
compromise between depth of field and focal point beam- 
width for each type of transducer. The transmit focal point 
in the elevation direction is controlled by a fixed-focus 
cylindrical acoustic lens. The lens is oriented such that it 
does not affect the azimuth focus. Hence, the elevation and 
azimuth focal points can be chosen to be different. 



Like the transmit focal point, the receive focal point in 
the azimuth direction is controlled by the delay given to 
each element of the phased airay. A basic difference is that 
because the echoes from Lncreaslng depths are recmved 
over a period of increasing time, the receive focal point 
can be controlled to track the depth from which the echoes 
are returning dviiainicaUy. This dynamic focttsing elimi- 
nates the depth of field problem associated with a fixed 
focus. 

The effect of dynamic focusing on image quality can he 
seen by comparing two images of a test target or "phantom" 
(Fig, 3). This phantom contains a series of nylon wires and 
cystic regions imbedded in a tissue-mimicking medium. 
The size of the wire images and the clarity of the round 
cysts provide a measure of azimuth resolution. In the first 
image the dynamic focusing is disabled and the receive 
focal point is fixed at 8 cm (Fig. 3a}. In the second image 
the dynamic focusing is enabled, providing optimum beam 
width along all the scan lines (Fig. 3b], 

The effect of elevation resolution on image quaht>^ is not 
as directly apparent as the effect of azimuth resolution. 
Elevation resolution determines the thickness of the sector 
image plane, AH echoes within this plane are mapped to- 
gether into the two-dimensional image. A thick image plane 
results in the obscuring of some small details. The image 
may still look sharp, but some information is lost. Consider 
an image of a small vessel in the body parallel to the plane 
of the scan. The cross section of the vessel should appear 
as a narrow black void because the blood within the vessel 
returns no echoes. If the image plane is thicker than the 
diameter of the vessel, then the echoes from tissue above 
and below the vessel axe mapped into the image of the 
vessel, obscuring it. This effect can be seen in two images 
of a tissue-mimicking phantom containing cylindrical re- 
gions representing horizontal blood vessels [Fig. 4)^ The 
images were made with two transducer.^ having different 
elevation focal points. 

Because there is no overall best choice of elevation focal 
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Fig. 2. Transmit beamwidth versus distance for equal-aper- 
ture, 5-MHz transducers wttti different focai lengths. 
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Fig. 3. Effect of dynamic focusmg 
on image quaiity. (a) Image using 
a transducer with a fixed recetve 
focus at 8 cm. (b) image using a 
transducer dynamicaily focused 
during reception of echoes. 



point, several transducers with different elevation focal 
points are available for ttie HP 77020 A Ultrasound Imaging 
System. (Table I). Thus, users can select the transducer 
frequency and focal point most appropriate for their par- 
ticular imaging application. 

Table I 

Transducers for HP 77020A Ultrasound Imaging System 



Transducer 


Center 


Elevation Focal 


Model 


Frequency 


Point 


21200A 


2.5 MHz 


12.0 cm 


21205A 


3.5 MHz 


7.5 cm 


21206A 


3.5 MHz 


4.0 cm 


21210A 


5,0 MHz 


6.0 cm 


21211A 


5.0 MHz 


3.5 cm 



Range or axial resolution is a measure of the ability to 
resolve individual closely spaced targets along the acoustic 
beam axis. Axial resolution in a pulsed echo system is 
determined by system bandwidth and center frequency. 
Transducers typically have a constant percentage band- 
width, and thus the higher the center frequency of the 
transducer, the wider the hand width and the better the 
range resolution. The bandwidth of the transducer is not 
the only consideration. The receiver bandwidth also affects 
resolution since the time response of the receiver filter can 
broaden the returning echo pulses. The 7 7020 A minimizes 
receiver pulse broadening by using a high-frequency IF 
stage and optimizing filters for each transducer center fre- 
quency. 



Side lobe Levels 

Another very important aspect of image quality is the 
effect of sidelobes and grating lobes. Sidelobe and grating 
lobe artifacts have plagued phased-array systems in the 
past. Sidelobes appear in an image as lateral streaks emanat- 
ing from ver^^ bright targets. The two images in Fig. 5 illus- 
trate the effect of sidelobes on image quality. The first image 
was made with sidelobes at a -40 dB level {Fig. 5 a], The 
second image was made with sidelobes at a - 60 dB level 
(Fig. 5b). The effect of the higher sidelobe level is the lateral 
smearing seen in the first image. Also, the dark blood ves- 
sels become filled with ^'clutter" caused by the higher side- 
lobes. 

The 7 7020 A suppresses the clutter level caused by high- 
order sidelobes by more than 60 dB and eliminates the 
grating lobe by setting the transducer element-to- element 
spacing to a half wavelength. 

Depth of Penetration 

One fundamental problem in ultrasound imaging is that 
soft tissue attenuates ultrasound roughly linearly with fre- 
quency. The human liver, for example, attenuates ul- 
trasound by about 0.5 dB/MHz/cm. Thus, a 2,5'MHz pulse 
propagating round-trip through 40 cm of liver tissue is 
attenuated by 50 dB and a 5,0-MH2 pulse is attenuated by 
100 dB. For this reason, deep imaging is usually done with 
lower-frequency transducers and shallow imaging is done 
with higher-frequency transducers to take advantage of 
their superior image quality. It will be shown later that a 
higher acoustic frequency does more for the image than 
just provide better range resolution. 

Frequency-dependent attenuation has an adverse effect 
on the spectrum of the transmitted pulse. Because high 
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Fig. 4. EffBCt of transducers eie- 
vation focat point on image quaiity. 
(a) image with an elevatfon focai 
point at 8 cm. (b) image with an 
elevation focal point at 4 cm. 
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Fig. 5. EffBCl ot transducer 
stdslobe fevBi on mage quafity. 

(a) image wiih sidelobe /eve/ at 
-40 dB Now the iateral smeanng. 

(b) Image with sideiobe /eve/ at 
-60 dB. 



frequencies are attenuated more than low frequencies, the 
center frequency of the transmitted spectrum is shifted 
down as the pulse propagates. Fig. 6 Ulustrates the shift in 
center frequency of a pulse spectruni at various depths. 
This shift in center frequency degrades the azimuth reso- 
lution because the wavelength of the pulse increases while 
the aperture size remains fixed. The effect of this spectral 
shift on image quality is controlled to some extent by the 
choice of receiver bandwidth. A wide bandwidth allows 
imaging at greater depths in the body, but with degraded 
resolution. A narrow bandwidth improves both azimuth 
and elevation resolution and focusing at the expense of 
axial resolution. 

The 77020A uses a compromise between these two ex- 
tremes for sector imaging, and uses a wider bandwidth for 
M-mode where axial resolution is more important. The 
ability to optimize filtering for each transducer and mode 
of operation is an important feature of the 7 7020 A System- 
Tissue Representation 

The ability of an ultrasound imaging system to represent 
the range of echo amplitudes as different gray levels in the 
image is another measure of image quality* This represen- 
tation affects how the separate tissues of the body are seen 



in the image. Echoes from within the body occur when the 
propagating ultrasound pulse encounters a change in the 
characteristic impedance of the body tissues. The charac- 
teristic impedance Z is defined as Z=pv, where p is the 
density of the medium of propagation and v is the velocity 
of propagation. When the impedance change occurs over 
a boundarj^ greater than a few wavelengths in extent, such 
as the boundary between a blood vessel wall and blood or 
organ tissue, a specular or mirror- like reflection occurs. 
The intensity I of echoes from these specular reflectors 
depends upon the reflection coefficient and is given for a 
boundary that is perpendicular to the direction of propaga- 
tion by I = ls(Z2 - Zi)/(Z2 + Z^) , where l^ is the incident pulse 
intensity, and Z^ and Z^ are the tw^o impedances. Specular 
echoes are frequency-independent and propagate in a di- 
rection that depends upon the orientation of the boundary 
v^ith respect to the pulse's propagation direction. 

Generally the impedance differences between the soft 
tissues of the body are small enough so that only a small 
portion of the incident pulse is reflected back to the receiver 
and the major portion of the pulse continues on to greater 
depths. However, there are cases where the impedance 
differences are great and most of the pulse is reflected back. 
When this occurs, a shadow appears in the image behind 
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Fig. 6. Variation with depth of the 
frequency spectrum of detected 
uftrasound echoes The incident 
pufse spectrum is centered at 3.5 
MHz (An attenuation of 7 dBl 
cm/ MHz is assumed.) 
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Scattered Size A 



Scallerer Size B 



Log Freqyency 

Fig. 7. Scatter intensity versus frequency for two scatterers 
of different size. The difference U-f^ between ttte scatter 
intensities at frequency f^ is greater than the difference /^-/j 
between f/ie scatter intensities at frequency /,. This result 
provides greater contrast between regions of tissue that are 
structuraHy different on a scafe much smaiier than the 
wavelength of the ultrasound putse. 

the impedance boundary. This effect can be most useful 
in applications where large impedance changes are abnor- 
mal, such as a calcification or gallstone in the gall bladder. 
The strong shadow behind a gallstone is often taken as 
evidence confirming its presence. 

When the region of an impedance change extends over 
less than one wavelength ^ part of the incident pulse is 
scattered in all directions according to Rayleigh scattering. 
The intensity of the scattered pulse is proportional to the 
product of the fourth power of the frequency of the pulse 
and the sixth power of the radius of the region of the imped- 
ance change. All soft tissues produce scatter echoes because 
of their cellular structure and small internal structural de- 
tails. Some of the scatter propagates back to the receiver 
and Is known as backscatter. When propagating through 
tissue, the transmitted pulse produces a continuous com- 
plex of backscatter echoes that constructively and destruc- 
tively interfere with one another in a random fashion. The 
resulting interference pattern seen in the image is known 
as a speckle pattern or simply as texture (see page 39). The 
texture intensity of each organ is different and can be used 
to differentiate organs and to indicate abnormal regions 
such as tumors within organs. 

The ability to show difference in speckle pattern inten- 
sities is another measure of image quality. Because of the 
frequency dependence of the scatter intensity^ higher-fre- 



quency transducers provide greater differentiation of tissue 
types. This is shown in Fig. 7 where scatter intensity is 
plotted against frequency for scatterers of two different 
radii. The improved tissue contrast detectability of higher 
frequencies can be seen in a series of three images of a 
tissue -mimic king phantom containing a region of scatterers 
made using three different frequencies (Fig. 8]. Notice the 
increase in detectability of the round 20-jum regions with 
an increase in frequency. 

Bandwidth also affects the appearance of the texture pat- 
tern. Since texture is caused by the addition and cancella- 
tion of randomly received echoes, a narrow bandwidth 
provides more coherence and produces a coarser texture 
than a wide bandwidth, The higher the frequency, the finer 
the texture, and the more filled-in or smoother the image 
looks, making it easier to identify boundaries between dif- 
ferent types of tissues or small vessels. 

Once the received echoes are optimally filtered, they 
must be detected without time-response degradation and 
presented as gray levels on a CRT display. The different 
types of echoes from the body cover a dynamic range on 
the order of 40 to 50 dB. That is, the weakest echoes can 
be 40 to 50 dB below the strongest echoes from areas of 
equal depth. Because the dynamic range of Images judged 
to be most pleasing to the eye is about 25 dB, a form of 
nonlinear processing is used to compress the range of inten* 
sities displayed in the image. This process is important to 
the quality of the image. Subtly different textures must be 
discernible while large specular echoes must not be too 
distracting. 

Two different methods for controlling the processing of 
echo amplitudes into image intensities are provided in the 
HP 77020A Ultrasound Imaging System's scanner. One 
method controls the blending of the logarithmic and linear 
versions of the echo amplitude before digitization and stor- 
age in memory. This variable blending allows the user to 
control the amount of echo dynamic range displayed in 
the image. For example, in abdominal imaging where soft 
tissue is to be displayed with high-level specular echoes 
from organ boundaries and bright gas-tissue interfaces* the 
blend is set to be mostly logarithmic. The other method 
digital !y controls the assignment to gray levels of the dig- 
itized echo amplitudes stored in memory. By using post- 
processing in series with logarithmic compression, the user 
can fine-tune the image presentation to provide the best 
differentiation between the tissue types studied. 




Rg. 8. Images made of a region of 20-fxm diameter scatterers surrounded by 7-ijun diameter 
scatterers using the 77020 A System's (a) 2.5'f^f-tz, (b) 3 5-MHz, and(c) 5.0-MHz transducers. 
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Coherent Speckle in Ultrasound Images 

by Paul A. Magnin 



SpecWe r3 the term used to describe the fnottling found in 
ultrasound images Specie nojse Tends to manites! itself as a 
relative ty high-contrasr, high-spat^aJ-frequency mask supenm- 
posed on the desjred ultrasound tmage Such nosse tends to 
make subtle drfferences in the gfay levels of different tissues 
imperceptible and reduces the apparent resolution of the image 
signfticantly below the diffraction limit. 

The speckle phenomenon is inherent ^n aJI coherent imaging 
modalities. Coherent illumination has the property that the 
phasors at ail points in the fieJd of illumination vary m unison 
While any two points in the field of illumination may have different 
relative phases, the phases of those points change in an identical 
manner The most commonty seen speckle is in laser images 
In this case the object illuminated with laser light has some sur- 
face roughness which is large compared to the wavelength of 
the laser light, The light returning to the eye of the observer from 
a particular spot on the object's surface is the compEex sum of 
the light returning from each infinitesimal surface of the tiny res- 
olution area. As long as the observer remains fixed relative to 
the object, the speckle pattern is stationary. If the observer 
moves, however, the speck Ee pattern changes- This motion re- 
sults in the scintillation the observer sees 

The reason speckle patterns are not observed in ordinary light 
is a result of the coherence time of the light, Sunfight and common 
artificial light are coherent, but only on certain time scales. That 
is to say. the relative phases of the illumination at any two points 
can vary in unison, but only for a very short time referred to as 
the coherence interval If one could make a photograph using 
ordinary sunlight in a single coherence interval, a speckle pattern 
would be obsen/ed The human eye. however, averages the 
swirltng speckle pattern to produce an apparently smooth image. 

Ultrasound speckle is caused by the constructive and destruc- 
tive interference of the backs catte red echoes. At any given Ume 
the echoes returning to the transducer emanate from a resolution 
volume. The resolution volume size is determined by the trans- 
ducer dimensions, the ultrasound pufse bandwidth and center 
frequency, and the degree of focus at the position of interest. 
The number of scattering sites within a resolution volume is typ- 
icaNy very large These scattering sites are normally individual 
tissue cells, tissue fibers, and collagenous connective tissue. 
Each scattering site in this volume produces an echo that can 
be thought of as a phasor which, when added to the echo phase r 
from all of the other scatterrng sites, produces the aggregate 
echo from that resolution volume. Att hough the echo from an 
Individual scattering sate is loo small to be delected at the fre- 
quencies typically used for medical imaging, the aggregrate 
echo from a collection of scattering sites can be detected If one 
assumes there are a large number of scattering sites in the res- 
olution volume and that ihe positions of these scattering srtes are 
randomly distributed wfth a uniform probability density function, 
then the amplitude and phase of the aggregrate echo from the 
resolution volume are random vanables The amplitude of the 
aggregate echo wiH be distributed according to a Rayleigh prob- 
ability density function and the phase of the echo will be dtstrtb- 
uied according to a uniform probability density function 

As the ultrasound pulse moves through the tissue, the number 
and orientation of the scattenng sites change, thereby changing 
the value of the complex? vector sum of the backscattered echoes 
The manifestation of this changing vector sum is the "salt and 




Fig. 1, Simp!e image (a) and compounded-ffequency imsge 
(b) of a tissue-mfmlckfng phantom made on a research 
phased-array scanner. 

pepper" appearance of ultrasound images. The random nature 
of the echo amplitude from tissue makes tissue identification 
from some form of echo signature very difficult, One cannot as- 
sume that a certain tissue will exhibit a characteristic echo. In- 
stead, some form of spatial averaging and statistical analysis or 
speckle reduction appears to be necessary before ilssue signa- 
ture analysis can be done 

A second vehicle responsibte for the mottled appearance of 
ultrasound images is the presence of inhomogeneous material 
whfch acts as an aberrating lens and results in some degree of 
phase cancellation This effect can be thought of as being analo- 
gous 10 looking through a translucent shower cunain. If one 
makes an image of a single point scatterer through an in- 
homogeneous medium, the phases of the echoes returning to 
each element of the transducer are disturbed from the spherical 
phase front one would othemise expect The resulting focused 
and detected echo has a randomly distnbuted amplitude 

One problem In medical imaging Is that some clinicians have 
a tendency to associate the speckle pattern with the tissue pathol- 



ocTOBEfl tsea Hewlett- packafid journal 38 



)Copr. 1949-1998 Hewlett-Packard Co. 



ogy. Although the speckle is a result ot scattering in the tfssue. 
the pattern itself is afso a function of tine imaging system. For 
example, the sjze of the "speckies" is related to the size of the 
resolution volume The resolution volume in turn is related to the 
transducer's aperture size, pulse frequency and bandwidth, the 
apod izat ion, and the nature of the subsequent signal processing, 
It is hard to imagine that any such pattem. which is so easily 
aftered by imaging system parameters, could be considered a 
function of the tissue alone. Optical texiure is commonly thought 
of as a surface variation with large spatial dimensions compared 
to the resolution of the imaging system. Similarfy. if an analogous 
ultrasound texture is to be examined, it is necessary to minimize 
the effects of speckle and to look at surface variations that are 
larger than the resolution volume of the imaging system. 

Two common approaches exist for reducing the amount of 
speckle artifact present in coherent images. The first involves 
spatially low-pass filtering Ihe image, which smoothes the mottled 
appearance of the image. However since the average speckle 
size is roughly equal to the resolution of the imaging system, the 
filtering also necessarily degrades the resolution of the system. 
For most applications this is considered loo high a pnce to pay 
The second approach involves varying certain imaging system 
parameters, which result in movement of the speckle pattern. 
and averaging the different speckled images Any parameter 
that alters the sum of the echoes returning from the tissue will 
after the speckle pattern. Phased arrays have an advantage over 
mechanically steered imaging systems in that the echoes re- 
ceived by the elements ot the array are summed to form the 
image Different amplitude weightings on the channels making 
up the array change the beam characteristics and alter the 
speckle pattern. One could also make images using different 
frequencies, foci, bandwidths, compression, or any combination 
of these parameters. If the different imaging conditrons selected 
produce independent or uncor related speckle patterns, upon 
summation the resultant composite image will exhibit a speckfe 
contrast reduction of 1/VN. where N is the number of images 
used to form the composite. Since fixed, resolvable structures 
remain in the same position for each ot the N images, their visibility 
relative to the speckle noise increases and no degradation of 
the system resolution occurs. Such images could be obtained 
in either a serial or a parallel fashion. The disadvantage of the 
serial implementation is the blurnng of rapidly moving structures 
such as heart valves. A parallel implementation eliminates blur- 
nng but requires increased system complexity and cost. 

One imaging system parameter which has been varied for the 
purpose of reducing speckle is the excitation burst center fre- 
quency^ In this scheme, referred to as frequency compounding, 
the transducer is excited with a gated sinu&oid, Each successive 



frame is made using a different sinusoidal burst excitation fre- 
quency. This is repeated for N different excitation frequencies 
and the IM video frames are summed or averaged visually by the 
eye of the observer Significant decreases In the speckle artifact 
can be obtained using this scheme. Fig. 1 shows two images of 
a tfssue-mimicktng phantom made on a research scanner. Fig. 
1 a is a conventional phased -array image showing a large cyst 
In the center and a smaller, barely visible cyst below and to the 
right of the large cyst. Rg 1 b is a frequency compounded image 
made using five different excitation frequencies. The speckle 
contrast is significantly reduced, leaving the cyst borders more 
clearly delineated. Furthermore, the small cyst is more clearly 
seen. Further reductions fn the speckle contrast using simitar 
techniques are necessary if more subtle differences in tissue 
texture and character are to be determined, 

Tissue CharacterfzatEon 

One potential diagnostic advantage of ultrasound imaging is 
its ability to differentiate between normal and pathological tissue. 
Presumably pathological tissue exhibits a different texture from 
normal tissue. Although one major difficulty m discerning subtle 
differences in tissue texture is the superimposed speckle artifact, 
it is possible that the speckle pattern itself may prove to be useful 
diagnostically. One could Imagine cases where tissue fibers are 
organized in some systematic fashion Speckle from such a tissue 
matrix would have a regular pattern and an orientation sensitivity. 
This is analogous to Bragg diffraction found in crystalline struc- 
tures However, to date, such tissue regularity has not been 
demonstrated. 

One of the most promising indicators of tissue character has 
been the sScpe of the attenuation coefficient versus frequency.^ 
Tissues tend to have attenuation characteristics that increase 
linearly with frequency. Recent studies suggest that an Infarcted 
myocardium has a significantly different slope ot attenuation ver- 
sus frequency than a healthy myocardium. Again it is necessary 
to account for or minimize the effects of speckle. Tissue charac- 
terization by ultrasound would provide a new dimension of diag- 
nostic information. Rather than simply obtaining structural infor- 
maiion, ultrasound has the potential for imaging pathology. Ulti- 
mately It may be possible to determine the location, severity and 
type of a tumor or cyst noninvasively, and thereby enhance the 
diagnostic utility of ultrasound Images. 
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